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CHAPTER I 


INTRODUCTION 

1.1 Deflected Turbulent Jets 

Almost all flow phenomena of practical engineering Interest are 
turbulent and do not lend themselves easily to experimental, analytical 
or numerical investigation. Deflected turbulent jets are no exception, 
involving fully three-dimensional flow phenomena. Power plant chimney 
plumes, cooling holes ejecting air from turbine blades, and the airflow 
around V/STOL aircraft in transition flight are examples of deflected 
turbulent jets. The present study emphasizes lateral injected jets into 
tubular cross-flow which may also possess swirl - a situation occurring 
in gas turbine combustion chambers, as seen in Figure 1, and a more 
complicated example of fully 3-D flowfields. When the disturbing jet 

is introduced into the crossflow, vorticity is added to the flow. The 
circulation of this vorticity stirs the two fluids together as the flow 
progresses downstream. Gas-turbine combustor designers utilize this 
mixing phenomenon to promote rapid combustion and evenly mix the 
products of combustion with cooler injected air, accomplishing this in a 
compact space. 1 


1.2 Combustor Flowfield Phenomena 

High- intensity combustion takes place in gas turbine combustion 
2 

chambers . Whereas most fuel burning takes place in the primary zone. 


the secondary zone Is where combustion Is almost completed, and the 
dilution zone is where further temperature reduction and uniformity is 
achieved. In both can and annular combustors, lateral jets of cooler 
air through round holes penetrate the flowfield. Some of these jets 
amalgamate with the swirl-induced central recirculation region of the 
primary zone and provide sufficient additional air for stoichiometric 
conditions. Other lateral jets provide additional air to help complete 
combustion in the secondary zone and to cool and evenly mix the products 
in the dilution zone in preparation for the introduction of the flow in 
the turbine. Similar problems arise in ramjet combustors. Clearly the 
turbulent reacting flowfield is fully three-dimensional; the combustion 
designer has a formidable problem in aerothermochemistry, and the 
research and development task is to provide a route which leads to the 
accomplishment of design objectives more quickly and less expensively 
than current practice permits. Recent conferences 2 *^ and textbooks 4 
extensively review progress and problems in this area. 

1.3 Present Study Objectives 

The focus of the present study is to characterize the time-mean and 
turbulent flowfield of a deflected confined turbulent jet, entering 
laterally into tubular cross-flow which may also possess swirl. Part of 
an ongoing research effort at Oklahoma State University, the ultimate 
aim is to provide a data base for turbulence model advances used by 
combustor designers. The specific objectives for this investigation 
included: 

1. Flow visualization was used to highlight the important 
features and structures of the deflected jet. Neutral ly- 
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buoyant helium-filled soap bubbles, smoke, and multi-sparks 
are to be used. Flow visualization was used to characterize 
one deflected jet without swirl (swirler removed) and with 
swirler angles of 45 and 70 degrees. 

2. A computerized data acquisition and reduction technique was 
developed and implemented. Because of the large number of 
data point locations and random nature of the fluctuating 
velocities, a high-speed analog-to-dlgltal data acquisition 
system with a computer-controlled probe drive was designed and 
constructed. 

3. A six-orientation single hot-wire technique was used to map 
fully the time-mean and turbulent velocity flowfield. The 
data acquisition and probe drive system was used to manipulate 
the probe and hot-wire voltages. Jet to cross-flow velocity 
ratios of 2, 4, and 6 were used with no swirl in the cross- 
flow. 

1.4 Outline of the Thesis 

The first chapter of this six-chapter thesis is the introduction. 
Deflected turbulent jets and their importance in engineering practice 
are briefly described, emphasizing the significance of deflected jets in 
combustor design. Finally the objectives of this study are stated and 
justified. 

Chapter II presents a historical perspective into research in 
deflected turbulent jets. The test section geometry and parameters 
investigated by other workers are outlined. Past work at Oklahoma State 
University into related combustor flowfield characterization is 
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discussed. 

A complete description of the experimental facility and measurement 
equipment Is presented In Chapter III. The data acquisition and probe 
drive system* des'gned specifically for this Investigation, Is fully 
discussed. 

The measurement techniques employed are discussed In Chapter IV. 
Flow visualization via bubble, smoke, and spark-gap techniques Is 
described. The chosen single-wire multi -position technique is described 
In light of previous workers^*** sensitivity analyses. 

Chapter V and VI are the concluding chapters of this effort. The 
former presents results of the flow visualization and measurement 
techniques. These results are discussed thoroughly; velocity plots are 
related to flow visualization photographs. Chapter VI emphasizes the 
main conclusions to be drawn from this Investigation. 

Tables III, IV and V present the time-mean velocities and 
turbulence quantities in tabular form. Figures 17-23 present flow 
visualization photographs; Figures 24-41 are two-dimensional plots of 
the time-mean and turbulent flowfield. 
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CHAPTER II 


BACKGROUND 

2.1 Review of Previous Studies In Deflected 
Turbulent Jets 

A number of experimental Investigations of the jet In a cross-flow 
have been reported In the literature. As early as 1948, Callaghan and 
Ruggeri 2 ^ examined a heated 200°C air jet directed normal to the wind- 
tunnel flow. The heated jet exhausted Into the confined tunnel through 
a sharp-edged orjflce. Jet velocities from approximately two to seven 
times the cross-stream velocity were Investigated for penetration depth 
and mixing. Thermocouple and pitot tube rakes were utilized well 
downstream of the jet inlet to record penetration depth and mixture 
temperatures . Ruggeri, Callaghan, and Bowden 2 ® extended this work In 
1950 to Include square and elliptical orifices. 

In 1952, Ruggeri 29 extended his contribution to include various 
angles (30, 45, 60, and 90 degrees) of Injection. The angles required 
that a tube be utilized to Inject the jet, as opposed to a sharp-edged 
orifice. Ruggeri used Schlleren flow visualization to confirm his 
measurements and commented on the wall effects of the wind tunnel. 

Jordlnson,®® 1958, was the first to determine experimentally the 
trajectory of the jet cross-section and defined the jet axis as the line 
connecting the points of maximum velocity. He also demonstrated that 
the cross-section of an Initially cylindrical jet Is distorted Into a 
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'horseshoe' shape by the cross-stream shearing action. 

Keffer and Baines 3 *, 1963, contributed some very carefully measured 
velocity data. More Importantly, they studied the turbulent structure 
In the deflected jet and showed that similarity for a reasonably small 
range of velocity ratios (4, 6, and 8) could be shown. Their work was 
the first use of a jet-oriented coordinate system, relating the jet 
centerline to that of a free jet. Pratt and Baines ‘ refined this work 
to account for scatter in the previous data, and also demonstrated that 
the profile of the jet Is conserved as a mixture between circular cross- 
section of random eddies and a pair of line vortices with small 
turbulences. Pratt aid Keffer 3 ^ continued this Investigation for 
various injection angles (60, 75, 90, 115, and 135 degrees) to the main 
flow. Finding differences between their jet trajectories and that of 
Keffer and Baines 3 *, they recognized the Importance of the jet Inlet 
velocity profile. In fact, their experiments used long tubular Inlets 
whereas Keffer et al used a simple orifice Inlet in the earlier study 3 *. 

Some of the first work to examine multiple deflected jets is 
represented by Norgren and Humenik 33 , 1968, intended to aid In the 
design of short-length combustors. As with the pioneering work of 
Callaghan^ 7 , they restricted their investigation to penetration depth 
and degree of mixing for heated jets. It should be noted that this work 
was one of the first basic research studies Into turbine Inlet 
temper ature profiles. 

In 1973, Campbell and Schetz 35 developed one of the first 
analytical models of a deflected jet and also verified this Integral 
method with experiments. Their model was successful In predicting flow 
velocities and trajectories in buoyant, heated jets. 
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Kamotonl and Greber 3 ® ware probably the first to study multiple jet 
Injection Into a confined cross-flow. The single row of jets was 
directed toward the opposite wall of a rectangular cross-section wind 
tunnel. Most Interesting of their conclusions was that the jets were 
only mildly affected by the opposite wall, unless the jet directly 
Impinged upon that wall. Holdeman and Walker 3 ^ extended this work and 
others to develop an empirical model predicting downstream temperature 
profiles. 

Chassalng, et al, 33 In 1974, contributed works comparing several 
zones of similarity of the jet In cross-flow problem. Krausch and 
Fearn* 0 contributed the first Investigation Into the properties of the 
vortex pair associated with the jet In cross-flow. 

The most detailed work thus far for a single deflected jet has been 
that of Crabb, Durao and Whltelaw**. Utilizing laser-doppler anemometry 
In the region of the jet, their measurements extensively quantify the 
velocity field with its associated turbulence and vortex pairs. The 
only detraction is that the jet to cross-flow velocity ratios of their 
measurements are quite low -- 1.15 and 2.30. Most investigators 3 ^" 33 
have been In the range of 4-10 times the cross-stream velocity as this 
approaches the dilution jet case more closely. The ratio of 1.15 would 
fall slightly above the region of film cooling. A tubular inlet was 
used to Inject the jet perpendicularly Into a large wind tunnel. 

The work of Rathgeber and Becker 39 Is representative of few 
Investigations dealing with jet injection into a cylindrical cross- 
section flow. They Investigated relatively small diameter jets as 
compared to the cross-flow diameter (cross-flow to jet diameter ratios 
between 17.2 and 50). These measurements quantify mixture and 
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trajectory data, but do not deal with turbulence details. 

Table I provides a ready comparison of the historical background of 
jets In crossflow. The table outlines parameters and techniques used by 
various Investigators and the variables measured. Note, for example, 
the number of Investigators using an orifice to inject the Jet where th» 
velocity profile would be hard to quantify. 

Research motivation Is provided from many sources. Claus*** points 
out that almost all previous Investigators except Crabb, et al 4 ~ have 
failed to fully report turbulence field. Inlet velocity profiles and 
vortex pair properties. His thesis Is that all of these data are 
extremely important to combustor modelers In confirming analytical 
tools. 

A recent (1979) review of the state-of-the-art in flowflelds 
modeling Is provided by LI 1 ley. 54 Holdeman and Srlnlvasan 3 present 
comments on NASA-Inspired modeling of dilution jets. References 49, 51, 
and 53 may give the reader some idea of the diversity of analytical 
flowfleld modeling schemes now under development. Additionally, 
Schetz^*, 1980, contributes a review of the entire field of Injection 
and mixing In turbulent flow. 

2.2 Past Work at Oklahoma State University 

Recently, as summarized by LI Hey 6 , experimental and theoretical 
research has been completed on 2-0 axlsymmetrlc geometries under low 
speed, nonreacting, turbulent swirling flow conditions, In the absence 
of any lateral jets. The flow enters the test section and proceeds Into 
a larger chamber (the expansion D/d * 2) via a sudden or gradual 
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expansion (side-wall angle $ » 90 and 45 degrees). A weak or strong 
nozzle may be positioned downstream to form a contraction exit to the 
test section. Inlet swirl vanes are adjustable to a a variety of vane 
angles with values of $ * 0, 38, 45, 60 and 70 degrees being 

emphasized. The objective was to determine the effect of these 
parameters on Isothermal flowfleld patterns, time-mean velocities and 
turbulence quantities, and to establish an Improved simulation In the 
form of a computer prediction code equipped with a suitable turbulence 
model . 

In conjunction with these research objectives, several experimental 
techniques have been developed Including: 

1. Flow visualization, achieved via still (for example, 
see Ferrell, et. al. 4 ^) and movie photography of 
neutral ly-buoyant helium-filled soap bubbles and smoke 
produced by an injector. 

2. Time-mean velocities have been measured by Yoon*** with a 
five-hole pitot probe at low swirl strengths. 

3. Turbulence measurements have been completed on swirling 
(up to $ * 70 degrees) as well as nonswirling flows using 
a six-orientation single-wire hot-wire technique by 
Jackson**, enabling all Reynolds stress components to be 
deduced. 

O 

4. An advanced computer code has been developed by Rhode 0 
and improved by Abujelala*^ to predict corresponding 
confined swirling flows to those studied experimentally. 
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5. Rhode’s® tentative predictions have now been supplemented 
by predictions made from realistic inlet conditions by 
Abujelala* 2 for a complete range of swirl strengths with 
downstream nozzle effects. 
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CHAPTER III 


EXPERIMENTAL FACILITY 

In many experimental efforts, ready-made components are either not 
available or prohibitively expensive. The facility described herein is 
a result of many man-hours of student time dedicated to design and 
construction. 


3.1 Wind Tunnel 

The test facility, dominated by the wind tunnel, is shown 
schematically in Figure 2 of Appendix 8. Air is induced to flow through 
a large foam inlet filter by an axial fan. The fan is connected by way 
of belts and pulleys to a seven-horsepower U.S. Motor varidrive. The 
vari drive permits the fan speed to be adjusted from 500 to 2850 RPM. 
The air is then forced into an expanding area section, where multiple 
20-mesh screens and a section of straws are employed to straighten the 
flow and significantly lower the turbulence intensity. 

Before introduction to the test section, the air flows through an 
ax i symmetric nozzle with an area reduction ratio of 25. The nozzle was 
built with a matching cubic radial profile, according to the method 
described by Morel 7 . The objective of this nozzle design is to produce 
a low turbulence level uniform velocity profile before introduction to 
the test section, with minimum adverse pressure gradient in the boundary 
layer. The exit throat diameter of the wind tunnel is approximately 
15 cm. 
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At the throat of the tunnel and before introduction to the test 
section, a variable-angle swirler may be fitted. The swirler consists 
of ten flat blades (with pitch-to-chord ratio 0.68) which may be 
individually adjusted to any angle from 0 to 70 degrees. Sander^ 
provides detailed information and measurements related to the swirler 
performance and swirler exit velocity profiles on this facility. 

3.2 Test Section and Dilution Jet 

The test section consists of a clear acrylic tube approximately 90 
cm in length attached to the wind tunnel throat. Standard commercial 
acrylic tube is used with 15.24 cm (6.0 in.) outside diameter, 0.318 cm 
(0.125 in.) wall thickness. The inside diameter is therefore 14.61 cm 
with a measured variation of ± 0.05 cm. To adapt the test section to 
the wind tunnel throat (inside diameter 15 cm), an adaptor section was 
machined to provide a smooth transition from wind tunnel thrc-at to test 
section. Two test section tubes were constructed. Both test sections 
have the dilution jet inlet located at x/D = 1.00 where x is measured 
from the tube inlet. The first tube has a series of probe access holes 
located at x/D = 1.00, 1.25, and 1.50 and at all six azimuthal locations 
270, 300, 330, 0, 30 and 60 degrees as shown in Figure 3. The second 
tube allows probe access to locations downstream of x/D s 1.50 (for 
example x/D « 1.75, 2.00, 2.50, 3.00) and at any azimuthal angle. This 
is accomplished via a tube rotation section, constructed from machined 
aluminum rings, acrylic, and ball bearings as seen in Figure 4. 

Laboratory compressed air at 6 to 7 atmospheres gauge pressure is 
used to supply the dilution jet. For stability, the supply air lines 


are large and are routed through two line regulators with an 
Intermediate tank (volume approximately 0.006 nr*) to dampen line 
oscillations. The second regulator was used to meter the flow rate. 
After the second regulator, the air was routed through a Fisher and 
Porter model 10A1735A rotometer for monitoring of the volume flow before 
introduction to the dilution jet. The dilution jet assembly, see Figure 
5, consists of a stagnation chamber, flow straightening section, and the 
jet nozzle. The stagnation chamber was constructed from 15 cm inside 
diameter aluminum pipe and filled with plastic household scrub pads to 
evenly distribute the internal flow. A hemispherical ly-shaped screen 
and convergent transition smooth the flow into the flow straightening 
section. Here the air flows through four brass screens for turbulence 
reduction. The nozzle was designed according to the method suggested by 
Morel 7 and is constructed of fiberglass. The nozzle diameter is 0.10 of 
the test-sect iun diameter for a cross-flow to jet area ratio (A c /Aj) = 
100. Construction of the nozzle, a multi-step process, consisted of: 
constructing a two-dimensional contour on a numerically-controlled 
milling machine from a computer-generated profile, using the contour and 
a hydraulic follower on a lathe to produce an axisymmetric male mold, 
and forming the fiberglass nozzle around the male mold, with aluminum 
flanges formed in. 

Once assembled, the dilution jet was attached to the air line 
coming from the rotometer and the nozzle was pressed into a special 
acrylic adaptor which is permanently attached to the test section. 
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3.3 Hot-Wire Instrumentation 


The sensing transducer used in this study is a normal hot-wire 
probe, DISA type 55P01. This probe has two prongs set 3 mm apart with a 
5 ym diameter tungsten wire between them. The exposed, effective length 
of the wire is approximately 1 mm, since the ends have been gold plated 
to strengthen the wire and reduce end effects. The probe support was a 
standard DISA 55H21 straight mounting tube. The anemometer used was a 
DISA type 55M01, constant-temperature standard bridge. The hot-wire 
voltage was measured with the computer-controlled data acquisition 
system, discussed in Section 3.5. 

3.4 Calibration Equipment 

A small axisymmetric free jet was employed to calibrate the hot- 
wire. The calibration jet facility consists of a contoured nozzle 
similar in shape to the dilution jet and wind tunnel contraction. A 
settling chamber and turbulence management section consisting of packed 
straws is just upstream of the nozzle. The nozzle exit diameter is 34 
mm. Using the standard laboratory air supply, the calibration jet is 
capable of producing Reynolds numbers up to 6 x 10 5 (based on throat 
diameter). The air supply is thermally stabilized by virtue of long 
indoor lines and is within i0.5°C of the facility temperature. The air 
is metered by means of a diaphragm valve and a Fisher and Porter model 
10A1735A rotometer. The jet was calibrated using a pitot probe 1.0 
nozzle diameters downstream of the exit plane. The temperature of the 
jet and the pressure before the rotameter were monitored during each 
calibration to account for minor variations from the initial calibration 
velocities. 
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The hot-wire was placed In the potential core of the jet during 
calibration. Utilizing a rotary table and two hot-wire support tubes 
(DISA 55H151, 01 SA 55H153), the hot-wire was calibrated in the u, v, and 
w directions as shown in Figure f. 

Jackson^ discusses in detail the merits of the chosen calibration 
expression 


E 2 = A + BZ 1/2 + CZ 

which is shown in Figure 7 for the three probe calibration directions. 
Figure 8 illustrates the pitch and yaw factors which will be discussed 
in Chapter IV. 


3.5 Data Acquisition and Probe Drive System 

The probe drive, shown schematically in Figures 9, 10, and 11, was 
specifically designed for these investigations. The probe is positioned 
in the flowfield by two stepper-motors, one motor for rotation and the 
other for translation. The probe is held within the square slider by a 
cylindrical holder with 0-rings to grasp the internal walls at any 
desired location. Both stepper motors step 200 times per motor shaft 
revolution. The rotation motor is geared down 3:1 so that 600 motor 
steps correspond to one probe revolution. For example, a probe rotation 
of 30 degrees requires exactly 50 steps. The software does not allow 
the probe to rotate In either direction more than one revolution to 
prevent cable twisting and coiling. 
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The translation motor is geared down 3:1 to a lead screw, which has 
a linear gear ratio of 2.24 revolutions per cm (equivalent to 5.69 revs, 
per Inch) translation. The effective step count for translation Is 
therefore 1344 steps per cm (3414 steps per Inch). With gear lash 
considered, the translation resolution is less than 0.03 mm and the 
rotation resolution is less than 0.5 degrees. The mass of the probe 
drive is approximately 3.9 kg (8.5 lbm) and is fastened to the test 
section with a large rubber binding strap. Reference 57 provides a 
detailed description of the probe drive and design philosophy. 

An Apple II computer was used to sample the hot-wire voltage and 
control the probe-drive stepper motors. A Burr-Brown SDM853 12-bit A/D 
converter was utilized to convert the 0-10 volt hot-wire signals to 12- 
bit digital words. The Apple II controls, via assembly code, the sample 
times and accepts the data as two 8-bit words directly in RAM. Further 
machine codes are used to reassemble the samples, take an average and 
standard deviation and store the results. BASIC code. Table IV, is used 
to reassemble these 8-bit words into decimal equivalent of 12-bit 
resolution. The system resolution is 2.44 millivolts. The data 
acquisition sample rate was fixed at 1000 samples per second for 5 
seconds. A higher sampling rate (up to 30 kHz) could be utilized with 
more memory available. Reference 56 provides a detailed explanation of 
the data acquisition system and the assembly codes used. 
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CHAPTER IV 


MEASUREMENT TECHNIQUES AND ANALYSIS 
4.1 Flow Visualization 

Flow visualization is used for primary Identification and 

characterization of the flowfleld, with three techniques being used. 

Bubbles, because of their reflective qualities and neutral buoyancy in 

the airflow, provide an excellent medium to determine the paths of the 

flow trajectories. Smoke, because of Its low comparative density and 

its tendency to mix well In the flow, makes an excellent medium to 
* 

follow the flow and to accent the turbulent paths and recirculation 
zones in the flow. A more novel flow visualization technique employed 
In these investigations Is the multi-spark method. With this technique, 
an Ionized path Is used to determine the relative velocity change from 
the position of the electrodes. The basic rig set up Is the same for 
all three flow visualization techniques, however differences occur in 
the type and quantity of lighting units, the camera time settings and 
light exposure times. 

4,1.1 Bubble Flow Visualization 

The bubble generator and Injection setup is shown In Figure 13, 
The bubble generator Is manufactured by Sage Action, Inc. It generates 
about 100 bubbles per second. The helium regulator range Is 0-207 kPa 
(0-30 psi) and flows at a maximum rate of 2.67 x 10" 6 m 3 /s. The maximum 
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bubble solution flow rate Is 2.50 x 10“ 8 m 3 /s and the maximum air flow 

rate Is 2.60 x 10“* mVs. A helium tank and an air line both with 

associated pressure gauges are connected to the bubble generator. The 

SAI bubble flow solution (BFS) Is Inserted directly Into a reservoir in 

the bubble generator Itself. The soap solution 1$ pumped out of the 

reservoir via helium. Three lines from the bubble generator are 

attached to the bubble Injector head. The head itself typically 

consists of three concentric tubes. The center one for helium, the 

middle annulus Is for soap and the outside annulus Is for air. Each 

line may be regulated by valves on the top of the bubble generator. A 

hole in the sleeve directly below the nozzle is where the injector head 

was inserted to inject bubbles parallel to the flow. 

•> 

The maximum bubble flow rate is 15.24 cnr/s. The slowest nozzle 
velocity Is about 4.2 m/s. From the equation Qy * Q n + where: 
where Qj « total volume flowrate, Q n * volume flowrate of the nozzle, 
and ■ volume flowrate of the bubbles. It was found that * 

0.0042; approximately one half of one percent. Hence the bubble 
flowrate has Insignificant effect on the nozzle flowrate. The Injector 
heads are also streamlined to minimize turbulent effects. 

Figure 13 shows the helium-filled soap bubble injection 
equipment. The lighting is approximately 3 m downstream of the test 
section. A light curtain dial provides light curtains from 0-1.5 cm 
wide through an adjustable slit and may be positioned to emit light 
angles from 0 to 360 degrees. The lighting Is on throughout the 
photography session and the exposure time Is determined by the camera 
settings. 
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The camera used was a Minolta SRI 200. The films used Include 
Kodak Trl-X Pan 400 ASA black and white, Ilford 400 ASA black and white, 
and Kodak color 1000 ASA film, with all of these giving excellent 
results. The camera was positioned approximately 0.5 m laterally from 
the test section and supported by a tripod. A low F-stop of 2 was used 
for maximum light Intake; the exposure time was set on B for a 5 second 
count. These settings were chosen after much testing to accentuate the 
bubble streakllnes Illustrating the flow trajectories. 

4.1.2 Smoke Flow Visualization 

The smoke generator and Injection setup Is shown In Figure 14. The 
generator Itself consists of a heating coil wrapped In steel wool and 
surrounded by a metal box with a drip tray In the bottom of the box to 
catch excess oil which may be drained out through a drain plug. 
Attached to the metal box Is a thermocouple which runs to a temperature 
Indicator outside the generator. Internal temperatures may then be 
monitored. The actual temperature of the heating coll may be 150°C 
greater than the temperature gauge reading. Experiments have shown that 
a temperature gauge reading between 250 and 300°C produces the optimum 
amount of smoke for flow visualization. It was also found that the 
temperature gauge must not exceed 350°C or a meltdown of the smoke 
generator gaskets will occur. The temperature Is adjusted through a 
rheostat which Is generally turned up to 50% power and then reduced to 
30% power for temperature stabilization. A valve above the generator 
may be opened to drip more oil onto the heating coll as needed for smoke 
generation. The air flow, which Is metered by the rotameter, runs 
through the smoke generator and up through the nozzle forcing the smoke 
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through the identical path. A few drops of oil produces a considerable 
amount of smoke and therefore has negligible effect on the previously 
monitored flowrate. 

The lighting device used for flow visualization consists of two 
commercial flash units. One unit Is a Vlvltar 2800 and the other a 
Sunpack 4220. The flash times may be adjusted on the units and the F- 
stop varied according to the distance of the camera from the flashes. 
Experimental results have shown that F-stop ■ 4 and flash time ■ 
1/2000 s with camera settings of F-stop equal 2 and exposure time from 
0.125 to 0.5 s produce the best pictures highlighting the flowfield 
features. The flash units are in a black box with a slit parallel to 
the test section and placed directly beneath and touching the test 
section. The positioning was chosen to minimize glare and maximize 
lighting through a vertical slit of light accenting a vertical cross- 
section of the flowfield. 

4,1,3 Spark-Gap Flow Visualization 

The spark-gap equipment schematic and configuration Is shown in 
Figures 16 and 17. The pulse generating circuit and pulse transformer 
Is manufactured by Sugawara Laboratories, Inc., Tokyo. The equipment Is 
capable of producing pulse trains of up to 200 pulses at frequencies 
from 1 kHz to 75 kHz. The output energy of the pulse Is 0.05 to 0.5 
Joules at voltages from 20 kV to 250 kV. As used in the present study, 
the electrodes are placed on opposite sides of the test section, 
typically one electrode above and one below the test section with a 15 
cm spark gap. Approximately 40 sparks are used with 0.5 J/spark at a 
voltage of 100 kV. Each spark pulse duration Is approximately 
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1 us ; time between pulses Is approximately 1 ms. 

When a high voltage source Is sparked across an air gap. an Ionized 
path Is created. Subsequent sparks will follow the current position of 
this low-resistance Ionized path. By judicious placement of the 
electrodes In the wall boundary layer, where there Is essentially zero 
velocity (next to wall), several discharges can follow the original 
Ionized path as It moves with the fluid. It Is necessary to have a low- 
conductivity test section material (for example, acrylic) so as not to 
Interface with spark paths. 

The spark Itself provides sufficient lighting for photographs. One 
camera (side view) Is used for photographs with zero swirl. Two cameras 
(side and end view) are used simultaneously In the swirl crossflow cases 
to give added perspective to the three-dimensional features of the 
technique. 


4.2 Quantitative Measurements 

In a turbulent, three-dimensional flowfleld the main flow direction 
may be unknown and conventional hot-wire or 2-0 Laser-Doppler techniques 
fall to supply sufficient velocity vector Information. To measure the 
three velocity components and their corresponding fluctuations, a three- 
wire hot-wire probe Is often used. Few 3-D Laser-Doppler systems are In 
use and are not cost-effective. 

As discussed by Jackson 11 , the three-wire probe technique has 
several drawbacks. Three anemometers are required. A multiple- 
orientation probe drive may be needed to align the probe with the mean 
flow direction. Because of the physical separation of the wires, 
spatial resolution of the probe Is poor. 
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Multi-orientation of a single hot-wire Is a novel way to measure 
the three components of a velocity vector and their fluctuating 
components. King*** modified a technique developed by Dvorak and 
Syred.*^ This method calls for a normal hot wire to be oriented through 
six different positions, each orientation separated by 30 deg from the 
adjacent one. Orientation 1 Is normal to the facility centerline, 
orientation 2 Is rotated -0 deg from this, etc. Mean *'nd rms voltcges 
are measured at each orientation. The data reduction 1$ performed using 
some assumptions regarding the statistical nature of turbulence, making 
It possible to solve for three time-mean velocities, the three turbulent 
normal stresses, and the three turbulent shear stresses. 

The six-orientation hot-wire technique requires a single, straight, 
hot wire to be calibrated for three different flow directions In order 
to determine the directional sensitivity of the probe. In the following 
relationships, tildes (~) signify components of the Instantaneous 
velocity vector In coordinates on the probe. Each of the three 
calibration curves Is obtained with zero velocity in the other two 
directions. The calibration curves of Figure 7 demonstrate that the hot 
wire Is most efficiently cooled when the flow Is In the direction of the 
u component (which Is normal to both the wire and the supports). The 
wire Is most Inefficiently cooled when the flow Is In the direction of 
the w component (which Is parallel to the wire). Each of the 
calibration curves follows a second-order, least-square fit of the form 

E 1 * A 1 + B 1 J 1 /2 * c i“i 
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which Is an extension of the familiar King's law. In this equation, A 1 , 
Bj, and are calibration constants and u^ can take on a value of 
u, v, and w for the three calibration curves, respectively. 

When the wire Is placed in a three-dimensional flowfleld, the 
effective cooling velocity experienced by the hot wire is 


">2 _ 2 ,w 2 
m + g Ir 


k 2 ; 2 


(4.2) 


where G and K are the pitch and yaw factors defined by Jorgensen®^ and 
deduced from the calibration curves. Those for this particular probe 
are given in Figure 8. Hence, equations for the effective cooling 
velocity can now be obtained for each of the six wire orientations. 
Simultaneously solving any three adjacent equations provides expressions 
for the instantaneous values of the three velocity components (u, v, and 
w in the facility x, r, and e coordinates, respectively) in terms of 
the equivalent cooling velocities. It is then possible to obtain the 
three time-mean velocity components and the six different components of 
the Reynolds stress tensor In the manner described by Janjua^ and 
Jackson. 11 

The uncertainty analysis Included a determination of the 
sensitivity of the six-orientation hot-wire data reduction to various 
input parameters that have major contributions in the response 
equations. Table II summarizes the sensitivity analysis performed on 
the data reduction program at a representative position in the swirling 
flow with $ * 38 deg. The table presents the percent change in the 
output quantities for a 1% change in each of the Important Input 
quantities individually, while the others are held at their standard 
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values. For the data presented in Table II, only quantities calculated 
from the probe orientations 1, 5 and 6 are used, for simplicity. This 
combination was chosen because the mean effective cooling velocity 
exhibited a minimum in orientation 6, and it Is expected 58 that in this 
case the combination 1, 5, and 6 will produce more accurate estimates of 
calculated turbulence quantities. The data of Table II demonstrate that 
the most serious inaccuracies in the measurement and data reduction 
technique are in estimates of turbulent shear stresses, the most 
inaccurate result being u V . 

Previously, in his measurements of strongly swirling vortex flows. 
King 58 compared his time-mean velocity and normal stress measurements 
with corresponding measurements obtained using a Laser Doppler 
Velocimeter (LDV). He found excellent agreement indicating the validity 
of the method. He was not able to compare shear stress measurements in 
his swirl flow, however, because he was unable to use his LDV for this 
purpose. In fact, despite the existence of advanced multicolor LDV 
systems, and their use for shear stress measurement, no one has yet 
reported their use in highly swirling flow situations: certainly not 

over a range of swirl strengths as reported in Reference 22. 

In the nonswirling confined jet case, results for time-mean 
velocities u and v, normal stresses u' r(T)S and v‘ rms and shear stress 
u'7 1 compare very favorably (see Reference 21, Figures 7 and 8) with 
those of Chaturvedi . 61 He used a crossed-wire probe for the shear 
stress measurements. So also did McKillop 14 for nonswirling confined 
flow In the same facility as Jackson and LI 1 ley. 22 Results, with and 
without exit nozzles, are In good agreement for the above quantities, as 
can be seen in Reference 14, Figures 21 through 28. 
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In the swirling confined jet case, 11 comparison with Janjua and 
McLaughlin®' 1 for a moderate swirl strength In an identical facility was 
made. They made triple-wire hot-wire measurements in a flow with an 
inlet swirl vane angle <t> - 38 deg., using analog-to-digital signal 
conversion and computer data reduction. For this purpose, it was 
necessary to know In advance the local time-mean velocity vector 
direction; the data of Yoon and Li 1 ley 1 ® (five-hole pitot probe) was 
used for this purpose. Their measurements® 3 of time-mean velocity 

compare very well with those of Reference 18 and hence of References 21, 
22, and 62. Measurements® 3 of the three normal Reynolds stresses and 
the three shear Reynolds stresses are compared at x/D * 0.5, 1.0 and 1.5 
with the six-orientation single-wire measurements of Reference 22. 
There is excellent agreement (see Reference 63, Figures 10 through 18), 
indicating again the validity of the present measurement technique. It 
appears to be an extremely viable, cost-effective technique for 
turbulent flows of unknown dominant direction. Probe interference 

appears not to be a major problem. Results are useful in recent 
prediction studies for confined swirling flows. 18 *® 4,65 

For the study of the technique presented by Jackson and Lilley 55 , 
Figure 5 through 9 of Reference 55 summarize measured values for the 
five representati ve situations in a turbulent flowfield. Each figure 
presents facility coordinate time-mean velocity, normal and shear stress 
values obtained with each of the five probe holder vs facility 

configuration possibilities of Cases 1 through 5. Case 1 is where the 
probe is in a nonswirling flowfield with the probe In facility 

coordinates. Cases 2 through 5 are where the probe is placed in 
different probe-to-facility orientations both in swirling and 
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nonswirling situations. A remarkable observation Is that, in general, 
the configuration Is of little importance, results appearing quite 
constant across the five cases. 

On the other hand, production run results^* 11 » 21 »22,62 have used 
the Case 1 configuration exclusively from each of the six possible 
combinations of three adjacent wire orientations. This was because of a 
lack of local flow directional knowledge; if this knowledge is available 

CO 

it is expected 30 that the combination with minimum cooling velocity in 
the central of the three wire orientations used will produce more 
accurate estimates of deduced flow quantities. In any case, the 
appropriate choice of wire orientation for minimum cooling velocity is 
not known a priori. However, for the turbulence quantities more 
confidence may be placed in the average of all possible wire 
combinations. This smoothing has been used exclusively in recent 
studies. 9 - 11 - 21 - 22 * 62 
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CHAPTER V 


RESULTS ANO DISCUSSION 
5.1 Flow Visualization 

Figure 17 shows very short time exposures of smoke tracing the 
extent of the deflected jet with jet-to-crossflow velocity ratio R of 2, 
4, and 6 in parts a through c respectl vely. The exposure time is of the 
order of 1 x 10"^ s, and vertical slit lighting is obtained with two 
commercial flash units. The camera is positioned to the side to obtain 
a view of the vertical rx-plane in the flowfield. Notice that as the 
deflected jet velocity increases, so does the jet penetration across the 
otherwise almost-paral lei crossflow. Clearly visible are the turbulent 
eddies - these are very structured near the injection location, and 
appear to extend further downstream in the jet direction at lower values 
of the injection velocity. Downstream of the deflected jet entry 
location. In the lower part of the main flow, a sequence of what appears 
to be eddies shed behind the lateral jet obstacle. A similar phenomenon 
was reported by Chassaing et al. 38 

Figure 18, 19 and 20 present long-time exposures of bubbles tracing 
the extent of the deflected jet with jet-to-crossflow velocity ratio R 
of 2, 4, and 6 in each figure respectively. Each figure presents swirl 
angles of 0, 45, and 70 degrees in parts a, b and c respectively. As 
with smoke tracing, the bubbles show the Increase of penetration with 
increase in R. The bubbles, however, show the time-mean boundaries of 
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the deflected jet. Clearly evident Is the lack of penetration In the R 
* 2 case; the R * 4 case crosses the test-section centerline at 
approximately x/D = 1,4 where the jet enters at x/D * 1.00 and then 
continues down the test section almost centered in the tube. The case 
of R 3 6 rapidly crosses the centerline (approximately x/D ■ 1.2, where 
the jet enters at x/D = 1.00), and continues downstream predominantly in 
the upper half of the test section. 

The swirl flow bubble pictures, shown in parts (b) and (c) of 
Figures 18, 19, and 20, illustrate the helical path of the jet and the 
strength of the precessing vortex core (PVC). In the case of $ = 70 
degrees swirl, sufficient negative axial velocity occurs in the PVC to 
carry the bubbles upstream to the swirler face, regardless of the R- 
value. For the cases of moderate swirl, $ = 45 degrees, there is a 
noticeable difference in the jet-to-crossflow interaction. The jet in 
the case of R 3 2 appears to mix broadly with the centrally-located PVC, 
indicated by the wide jet outline and broad PVC. The cases of R = 4 and 
6 exhibit less Immediate mixing with the precessing vortex core, tending 
to disturb its presence as the lateral ly-injected jet penetrates across 
the central part of the main flow, where the PVC would otherwise be. 
The R = 6 case in particular seems to "wrap around" the central axis - a 
smearing of bubbles can be seen on the tube inside wall. For the case 
of strong swirl, 0 = 70 degrees, there is very little difference in the 

flowfield between the different injection velocity cases. The cases of 
R « 4 and 6, however, can be seen to slightly deflect the swirl axis at 
approximately x/D ■ 2.00 (jet enters at x/D * 1.00). The case of R * 6 
does seem to lower the upstream penetration of the PVC, as exhibited by 
a lack of bubbles. 
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Figures 21, 22, and 23 present spark-gap flow visualization 
pictures for the same cases of R - 2, 4, and 6, using the method 
described In Chapter IV. These particular photographs were taken with 
the electrodes positioned at x/D * 1.50 where the jet enters at x/D = 
1 . 00 . 

In part (a) of these figures, the camera Is positioned to the side 
of the facility and a vertical rx-plane Is observed. In the swirl flow 
cases of parts (b) and (c), a second camera was simultaneously operated 
from a downstream location to illustrate the e-plane behavior of the 
sparks. In these swirl cases, both photographs have been combined to 
form a common picture. The respective cases R * 2, 4, and 6 with no 
swirl exhibit the change in the flowfield from x/D = 1.50 and continuing 
downstream. The case of R * 2 shows how the flowfield is merely 
deflected upward by the entering jet. The lower part of the arcs 
apparently are deflected around the jet, away from the control plane, 
hence a true 3-D effect on the photograph. The case of R = 4 shows 
flowfield acceleration above and around the jet, which has Its 
centerline nearly corresponding with the centerline of the tube. The 

'fold-over' just above the jet centerline probably corresponds to the 
downflow around the jet as the jet displaces the crossflow in the upper 
half of the test section. The case of R ■ 6 shows less uniform 
behavior. The arcs do appear to define the upper bounds of the jet and 
the turbulent region behind the jet. 

The swirl flows presented in Figures 21, 22, and 23 parts (b) and 
(c) are actually two photographs taken simultaneously by separate 
cameras; the two negatives are combined to print a common picture. 
Again, the electrodes were placed at x/D * 1.50, the lateral jet 
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entering at x/D = 1.00. A wire was placed in the centerline of the tube 
to prevent the spark from arcing to the tube walls and to help define 
the tube centerline. The end views exhibit a great deal of reflection 
off of the inside acrylic tube walls. 

With moderate swirl ($ * 45 degrees) the cases of R * 2 and 4 have 
little affect on the swirl pattern shown with this technique. The swirl 
pattern in the case of R ■ 6, however, is seen to be deflected by the 
jet. The swirl strength seems to be enhanced in the lower part of the 
test section by the additional momentum of the deflected jet. With 
strong swirl ($ =70 degrees), the cases of R = 2 and 4 seem to 
slightly inhibit the swirl strength, whereas the R * 6 case appears to 
have little effect except to organize the swirl pattern. 

5.2 Hot-Wire Measurements 

The time-mean velocity and turbulence quantities for the three jet- 
to-crossflow velocity ratios R * 2, 4, and 6 are presented in Figures 24 
through 41. The situation with jet-to-crossflow velocity ratio R * 2 is 
shown in Figure 24 through 29. Figure 24 has traverse angle 8 = 270 
degrees, Figure 25 has traverse angle 8 = 300 degrees, etc. Figures 30 
through 35 are for jet-to-crossflow velocity ratio R * 4; Figures 36-41 
represent ratio R * 6. Each figure is composed of twelve plots, (a) 
through (1), of the data for one traverse angle e at all seven axial 
locations (x/D = 1.00, 1.25, 1.50, 1.75, 2.00, 2.50, and 3.00). 

Using Figure 30 as an example, it can be seen that the time-mean 
and turbulent flowfields for the case of jet-to-crossflow velocity ratio 
R = 4 and traverse angle 8 = 270 degrees are presented. Recall that 
Figure 3 provides the geometrical relationship between the jet and the 
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traverse angle e. For e * 270 degrees, as In Figure 30, the viewer Is 
seeing an rx-plane of the flowfleld which passes through the test 
section centerline and Is normal to the lateral jet nozzle centerline. 
The top (bottom) of each plot corresponds to the first (last) measuring 
station, as shown In Figure 10. Tables III, IV, and V present the 
actual numbers used to produce the plots. 

Subparts a, b, and c In each of the Figures 24 through 41 present 
the normalized time-mean velocity component magnitudes u7u Q , 7/u Q , 
and w/u o respectively. Subparts d, e, and f give the normalized 
fluctuating velocity components (u^ ms /u Q , v^ s /u Q , and 
multiplied by 2. Subparts g, h, and i exhibit the three shear stresses 
(u'v 1 , u V , and v V ) normalized by u Q ^ and multiplied by 2 for 
plotting. Subparts j, k, and 1 provide the total velocity Y * 

o o o 1 / ? 

u + v + w ) ' , the axial turbulence Intensity normalized by the 

local mean velocity u^/u and the normalized turbulent kinetic energy 

A (ulL ♦ + w iL)/ u n 2 respectively. 

2 rms rms rms' o 

The plots were produced on a Tektronix 4006 terminal connected to 
an IBM 30810 using PLOT 10 as the graphics control language. The data 
are merely scaled and plotted point-to-point for each axial location. 
The x/D scales also provide as the magnitude scale for each normalized 
data point. For example. In Figure 24a, the values of Tf/u 0 at x/D = 
1.00 are scaled such that values of ZT/u 0 ■ 1.00 are placed at x/D = 
1.25. In this figure, the values of u/u 0 are very nearly 1.0 across the 
traverse except near the centerline, where the flow decelerates just 
because of the jet, which Is slightly below this traverse for R * 2.0. 
In Figure 24b, the values of V/u Q are much less than 1.0, approximately 
0.15, and are plotted as such. 
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5.2.1 Jet-to-Crossflow Velocity Ratio R ■ 2.0 


The velocity ratio R * 2 Is represented In Figures 24 through 29. 
Figure 27 provides the best perspective to visualize the flowfleld: the 
rx-plane shown (0-0 degrees) In It is the same as used for the flow 
visualization photographs. Figures 17. 18. and 21. Figure 27a shows how 
the jet, which enters from the bottom of the plot at x/D ■ 1.00, affects 
the axial velocity profiles. Figure 27b, the radial velocity plot. Is 
Interesting In that the six-position technique 4 s capable of accurately 
measuring the jet velocity as It Impinges on the probe parallel to the 
probe axis. Note that the techniques can only produce magnitudes; there 
are no negative 7 velocities for instance. In Figure 27a It can be seen 
that the R * 2 case has virtually no effect on the mean flowfleld above 
the centerline; this is borne out by Figure 24a, the axial velocity 
flowfleld in the horizontal rx-plane across the centerline (traverse 
angle 6 * 270 degrees). 

The normal stresses, see parts d, e, and f of Figure 27, are spread 
through the flowfleld over a wider region than the mean velocity. These 
and the shear stresses in parts g, h, and 1 are very low In magnitude. 
Part 1 shows that the region of significant turbulent kinetic energy 
extends only a short distance. The total velocity magnitude, Figure 
27j, Is well mixed and evenly distributed across the test section by x/D 
* 3.0. These observations for the R ■ 2 situation are not appropriate 
at higher values of R. 

The traverse angles e 3 330 degrees and e ■ 30 degrees, show no 
particular surprises for the R * 2 case. Figure 26a, the axial velocity 
In the rx-plane 30 degrees from the vertical shows an Interesting 
acceleration and then (at x/D » 2.0) deceleration In the mean flow 
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affected by the lateral jet. Figure 28a shc-'s an almost Identical 
configuration, Indicating good symmetry about the vertical plane. 

5.2.2 Jet-to-Crossflow Velocity Ratio R ■ 4.0 

The jet-to-cros$flow velocity ratio R * 4, as can be seen In Figure 
19, provides a flowfleld which Is more Intricate than was the case when 
R ■ 2. The centerline of the jet crosses the centerline of the test 
section, 43 and smoke flow visualization photographs, given In Figure 17, 
show that the turbulent eddies are large compared to those of R * 2.0. 

Focusing attention on Figure 33a, which shows profiles in the 
vertical plane through the centerlines of both jets. It can be seen that 
the jet has a marked effect on the mean flowfleld below the jet 
centerline, causing axial flow deceleration. The axial velocity at x/D 
* 1.00 shows that the upstream flow has slowed to go around the lateral 
jet on either side; the velocities on either side are greater, as seen 
In Figure 35a (e * 60 degrees) or Figure 31a (e » 300 degrees). 

Because the jet centerline crosses through the facility centerline, 
Figure 30 (8 * 270 degrees) provides insight into the structure of the 
jet. First, parts a, b, and c show good symmetry. The jet centerline 
was previously measured to cross the crossflow centerline at x/D * 1.35 
by Ferrell, et al . 43 This is seen to affect the tangential and radial 
velocities first at x/D ■ 1.25. The axial velocity Is not changed until 
x/D « 1.50 where it exhibits a flat acceleration and then, at x/D « 
1.75, a deceleration in the mean velocity. This deceleration 
corresponds to an increase in the surrounding axial velocity outside of 
the jet. By x/D ■ 3.00, the axial velocity profile is relatively flat 
again. In reviewing the tangential velocity profiles in part c of the 
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figure, the most surprising feature Is the symmetry and uniformity of 
the curves. At x/D * 1.75, 2.00, 2.50, and 3.00 the "gull -wing" shape 
of the profiles Is very likely caused by the data reduction results. 
That Is, the profiles should actually look like "normal probability 
distribution curves" — If negative values could be obtained. 
Physically this means that the tangential velocity along the central 
plane of the jet first rotates one direction, say clockwise, outside of 
the jet, and then the other direction, say counter-clockwise, Inside the 
jet. The normal stresses. Figures 30 parts (d) and (e), exhibit the 
expected peaks at the crossover. Note the dual peak In the u-dlrectlon 
normal stress. Figure 30d, at x/D ■ 1.75. 

There are some problems with the measurement technique and 
assumptions In signal Interpretations In some regions of the flow. The 
normalized axial velocity In Figure 33a, shows a large value below the 
jet centerline, x/D ■ 1.25. The shear stress u , w , /u Q shows a very large 
and erroneous value at the same location, and is related to the 
erroneous value of the mean axial velocity. These erroneous values have 
been faithfully presented along with the rest of the data. 

5.2.3 Jet-to-Crossflow Velocity Ratio R « 6.0 

The jet-to-crossflow velocity ratio R » 6 case dramatically 
exhibits some of the behavior expected of these deflected jets. 
Examining Figure 36a, which gives data in the horizontal plane, It can 
be seen that the jet centerline crosses the facility centerline at 
approximately x/0 » 1.15 to 1.25, The axial velocity has a marked 
depression inside the jet and a large acceleration around the jet sides, 
as If the crossflow were passing a solid body. Figure 36c Is quite 


Interesting in that it appears that all of the profiles could cross the 
zero line If the method would account for negative values. 

As for the R * 2 and R • 4 cases, the plot that displays the data 
next in an interesting format Is the traverse in the vertical plane, 
with e*0 degrees, Figure 39. As with the previous two cases, the 
axial velocity profile can be used to locate the jet centerline. 
However, the total velocity profile, given In part (j) of the figure, Is 
actually more accurate In locating the maximum velocity centerline. 
Providing testimony to the accuracy of the technique Is Figure 39b, 
where the radial velocity 7/u Q Is seen to asymptotically approach the 
normalized lateral jet Inlet velocity as the probe Is lowered toward the 
jet exit throat. Unfortunately, the shear stress plots, given In parts 
(h) and (1) of the figure, show that the technique is very sensitive to 
erroneous readings and probably to dwell time. The sensitivity of this 
technique to Input variables was discussed In Chapter IV and was the 
subject of analysis by Jackson. ^ 

5.2.4 Assessment of the Measurement 

The six-position hot-wire technique is remarkablj reliable In 
measuring the the time-mean velocities for the flowfields In this 
Investigation. As evidenced In Figure 24, 30, and 36 for 6 * 270 
degrees, It can be seen that the technique is repeatable In terms of the 
flowfleld symmetry. These figures display the lack of probe 
Interference effects as well — the same measurements are obtained on 
either side of the symmetry plane. Most surprising Is the fact that the 
technique can resolve the component direction even when that direction 
Is normal to the wire In all six orientations as also found in a 
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directional sensitivity study.® 5 The normal stresses are reasonably 
reliable In that there are few large discontinuities In the data. The 
shear stresses* however, exhibit less continuity, but the 
discontinuities do exist In regions of large shear such as behind the 
jet. 

By using the maximum velocity magnitude to define the jet 
centerline, Figure 42 displays the comparative centerline locations for 
these experiments as compared to the Infinite crossflow situation.^® As 
expected, the jet penetration for the confined cylindrical situations Is 
reduced from that of the Infinite crossflow situation, although low 
values correspond to jet injection velocities for which the confining 
boundaries have little effect. 
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CHAPTER VI 


CLOSURE 

6.1 Conclusions 

Experiments have been conducted to characterize the time-mean and 
turbulent flowfield of a deflected turbulent jet in a confining 
cylindrical crossflow. Jet-to-crossflow velocity ratios of 2, 4, and 6 
were investigated, under crossflow inlet swirler vane angles of 0 
(swirler removed), 45 and 70 degrees. Smoke, neutral ly-buoyant helium- 
filled soap bubhjes, and multi-spark flow visualization were employed to 
highlight interesting features of the deflected jet, as well as the 
trajectory and spread pattern of the jet. A six-position single hot- 
wire technique was used to measure the velocities and turbulent stresses 
in the nonswirling crossflow case, as a demonstration of improved data- 
acquisition capability. A computerized high-speed data acquisition and 
probe drive were designed and constructed to manipulate the hot-wire and 
reduce the varying voltages to the statistical mean and root-mean-square 
voltage. The voltages were then reduced to the time-mean velocity and 
turbulent Reynolds stresses with a Fortran computer code. 

The high-speed data acquisition system enabled three entire 
flowfields to be characterized for time-mean velocities, normal and 
shear stresses, for three different lateral jet injection velocities 
into nonswirling crossflow. The multi-orientation technique worked well 
for time-mean velocities, normal stresses and most of the shear 


37 


stresses. The extensive results are printed in tabular form and 
presented in rx-plane plots, in a manner useful to flowfield modelers. 
As expected, measurements confirmed that the deflected jet is 
symmetrical about the vertical plane passing through the crossflow 
axis. The jet penetration into the nonswirling crossflow was found to 
be reduced from that of comparable velocity ratio infinite crossflow 
cases. The flow visualization techniques enabled gross flowfield 
characterization to be obtained for a range of lateral jet-to-crossflow 
velocity ratios and a range of inlet swirl strengths in the main flow. 
The swirl in the confined crossflow was found to deflect the jet from 
its vertical course in a spiral fashion. However, the jet still gets 
absorbed finally into the precessing vortex core (PVC) of the 
crossflow. Evidence was also found that the jet can deflect the axis of 
the PVC and hinder the upstream propagation of it toward the swirler. 

6.2 Recommendations for Further Work 

The stage is now set for even more complete surveys of the time- 
mean and turbulence properties of deflected jets. Fundamental research 
should be continued first with the addition of a second jet directly 
opposing the original, and later with multi-jets at one axial station, 

1 O 

to complement NASA-Lewis work A » J in rectangular ducts. For ease of 
representation, a fully three-dimensional plotting technique should be 
developed and implemented in addition to streamlining of the data 
reduction technique. Equally important to more complete flowfield 
investigation will be the capability to measure the turbulent 
dissipation rate. A computer with larger RAM (random- access memory) and 
faster clock speed would enable rapid signal analysis and interactive 
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dwell time estimation to reduce further measurement time, as well as 
enable dissipation rate measurements. 
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APPENDIX A 


TABLES 
TABLE I 

PREVIOUS INVESTIGATIONS OF JETS IN CROSS-FLOW 


AEf. 

JET 

diameter 

(an) 

INCIDENT 

ANCLE 

JET 

inlet 

STYLE 

CROSS-FLOW 

VELOCITY 

(*/») 

VELOCITY 

RATIO 

GEOMETRY 

MEASURED 

TECHNIQUE 

27. 

6.35.9.5 

12.7.15.9 

90 

orifice 

- 

2-7 

open 

penetration 

parameters 

T/C and pitot 
rakes 

26. 

— 

90 

square, ellip- 
tical orifice 

46 

2-6 

open 

penetration 
and mining 

T/C and pitot 
rakes 

29. 

6.35.9.5 

12.7,15.9 

90,60, 

45,30 

pipe 

71.6,121.9 

2.9-5. 7 

open 

penetration 
and mixing 

T/C and pitot 
rakes 

3C. 

12.7,25.4 

90 

orifice 

- 

4.6.8 

open 

total press., 
flow direction 


31. 

9.5 

90 

pipe 

1.5 

4.6.8 

open 

velocity, turb- 
ulence Intentlty, 
entrainment 

oriented hot- 
wire 

32. 

6.35, 9.5, 
12.7 

90 

orifice 

0.914,3.66 

5,15.25.35 

open 

profiles and 
penetration 

photographs 

33. 

— 

90 

pipe 

-* 

0.55-2.20 

multiple 

Jet* 

penetration 
tnd mixing 

T/C and total 
press, probe 

34. 

6.35 45.60.90. 

105,120,135 

pipe 

1.56 

4.6.8 

open 

trajectories, 

velocities 

hot-wire 

3S. 



ANALYTICAL MODELING, NO EXPERIMENTS 



36. 

6.35 

90 

pipe/nozxle 

6-9 

2. 6-8. 5 

tquare no let 
confined 

velocity and 
temp, profiles 

T/C and hot- 
wire 

37. 

6.35-26.4 

90 

orifice 

15 

1.67-5.67 

2-dlmentional 

penetration, 

mixing 

T/C rakes 

38. 

40.0 
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pipe 

3.4 

2.37,3.91, 

6.36 

multiple 

Jett 

velocity, temp. 

similarity 

profiles 

hot-wire 

39. 

3.23.4.57 

6.30.9.32 
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pipe 

6.15 

2.4-12.4 

confined 

cylin- 

drical 

channel 

penetration, 
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101.6 
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90,105 

pipe 
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4.8 

open 

vortex 

strength 

pitot 

probes 

41. 

25.4 

90 

pipe 

12 

1.15.2.M 

open 

detailed velocity LDA, X-wire, 


and turbulence heliun trace 
field 


42. 


ANALYTICAL MODELING. NO EXPERIMENTS 
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TABLE II 


EFFECT OF INPUT PARAMETERS ON TURBULENCE QUANTITIES IN THE 
SWIRLING FLOW WITH $ - 38 DEG. AT A REPRESENTATIVE 
FLOWFIELD POSITION (x/D = 1, r/D « 0.25) 





% CHANGE IN 

TIME-MEAN AND TURBULENCE QUANTITIES 



t change in 










PARAMETER 

PARAMETER 
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r ms 

w' 

rtns 
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vV 

T 1 
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? 5 
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♦12.88 
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♦14.29 

♦ 7.46 

? 6 
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-10.59 

-0.36 

-8.50 
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♦7.07 

-9.54 
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-54.29 

-11.94 

E U 

♦ 1 
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-0.06 

♦0.14 

♦ 1.63 
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♦ 0.39 

fl 

♦ 2.86 

♦ 1.49 


♦1 

♦0.05 
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♦ 0.14 
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-0.13 

♦1.57 
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♦ 1.49 

E L., 
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-0.16 

♦ 0.18 
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-0.63 

♦ 1.03 

-1.08 

-2.0 

-5.71 

0.0 

G 

♦1 

-1.02 
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-1.01 

-1.0 

0.0 

-0.98 

-2.0. 

-2.86 

-1.49 

K 

♦ 1 

♦0.01 

-0.04 

♦ 0.01 

♦ 0.01 

0.0 

♦ 0.01 


O.C 

0.0 

% z 0 

♦ 1 

♦ 0.05 


♦ 0.14 

-0.13 

-0.13 

-1.77 


Ai) 

CD 

Csl 

♦ 1.49 

Vr 
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♦ 0.21 
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♦0.13 

1 

o 

vO 
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Tz p z R 
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-0.16 

1 

-0.08 

♦0.13 

0.0 

4-0.69 
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-2.0 

♦2.86 
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TIME-MEAN AND TURBULENCE DATA FOR JET TO CROSSFLOW 
VELOCITY RATIO R = 2.0 
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TABLE III (Continued) 
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TIME-MEAN AND TURBULENCE DATA FOR JET TO CROSSFLOW 
VELOCITY RATIO R = 4.0 
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TABLE IV (Continued) 
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T ABLE IV (Continued) 
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TIME-MEAN AND TURBULENCE DATA FOR JET TO CROSSFLOW 
VELOCITY RATIO R = 6.0 


m n m 

oo in o 

0 ) n in m h- 

00 

CD CO CD 


h* 

on ix> in O I- 

in oo in n 0 ) 

* 04 

O * 3> 


CO 

h- O 0) 

o 

04 CO C7) CD 

in ^ in (j) c* oi 

0) n O 


04 CO ID 

00 

u» <n p 

CN 

CD 

Omonri 

» M B ' 10 O 

in ix> O 

MT O) CD 

r** 0) O m ^ ^ 

CO 

r* O 


0) 


0»-'Q¥® 

000600 

O O - 

'OO 

OUlOOOO 

o 

' 0) 0) o 

oo 

o 

0) 0) 0) o O 


^ o - - - - -oo-" -ooo 



^ 04 a) in 

r> o CD n o oo 

0) co O t£ n a 

o) 04 m in co a) 

— CM 0) CM V O 


ninBinnO 

n co co in 

O* CO 0) (7) CO O 

cm O n 

— r~ co — r- cl 

o 

in io n O 0) n 

r* cn m co 

in ’t O in in o 

n S- o w r> n 

OIMO-CMI- 

in 

oo u> 0) in 04 o) 

O o io O — io 

r B O H h* n 

tt a) ro 

0* ID 04 CD 

04 

000"0 

- O 0> O O O 

<7> <7> <7) (7) O 

— 0> 0) 0) 0) o 

o o ai at o o 


-’-O’-’-- 

-oooo- 

-oooo- 

- - o o - - 


04 O 04 00 00 

io cm in o io 0) 

04 O* ^ O* CD ^ 

0) CM O — 1- CM 

r* 04 CD CD in 04 


0) 0) o oo co in 

oo — r- 0 ) in a) 

in — oo in 0) 

oo ^ o- cc oo in 

CO ^ h- u> V ^ 

Q 

O io f- — o n 

O CD CD O <7) <7> 

* ID 0) - 0) - 

nnip^NO 

m co 0) to 0) 

O 

r- o o> io O to 

O <7 h- 04 0 ) 0 ) 

— O co co © 

04 CO n in 04 TJ 

CO CO in (7) <7) *7 

04 

O - — cm n - 

o - o o - - 

— — 0) CO O — 

- O eo f~ ai - 

O O ® f- « o 



- - OO - - 

——OOO— 

- - o o o - 


— w in — in — 

O O CO * 10 - 

04 in in a) n 

(0 04 tf) (7) co ^ 

O ri a cm b in 


00 00 0) 04 O) 0) 

co It) O V) - 0) 

in ^ o) 0 ) 04 co 

in <7 04 0) 04 

u> n o ® n m 

in 

eo in o- 0 ) — 04 

O ® in i- u cm 

— oo ^ in m a) 

f- CO rr in TT (7) 

- ® S ® a> v 

o* 

O *"'*"*" n n 

CD C7) 0) in CO •“ 

in in oo a) 04 in 

CO O c- CO CD — 

rt in co in oo 

• 

04 04 

--BBBN 

— CM — O CM CM 

04 0) 00 O) 04 

co o- co O 





——OOO— 

——OOO— 

\ 

X 

n OO 04 04 in 

m o* n co <r> co 

CO ^ (7) CD ID C7) 

r» in h- in O id 

o> ip m ^ ^ o* 


in eo ■<* — (*) o 

10 9 O 0) 10 10 

OcooioiOn 

oo in r- in ^ ^ 

O O n n — ® 

o 

in oo in n ^ 04 

04 04 tt co co co 

co in o 0 ) co in 

r- 04 r** — in r- 

-oinmapi 

in 

n oo ip id n h 

in * a) 04 in 

O co in in cm — 

aiONMOin 

® U> <T) 0) O O 


O 0 ) W 0 ) o 

- - o - o - 

04 CM CD ^ CD 04 

04 n n in in co 

CD 04 CD *7 


"OOO' 






m *- co h* r- 0 ) 

in 0 ) co o i- i” 

CO <7 (7) in 04 (0 

0 ) m 04 m to in 

OM»0)CMN 


oinn^j m 

cm c- 0 ) r~ O in 

— u) o ® in co 

in in m 04 <7) 

CM n ^ o CM ® 

in 

. . * in rr in *7 

<D O ’J f) r~ — 

tt co 04 in — in 

in co o in 04 

cobnixk 

04 

— CM O t- 0) CM 

n n <7) in co 04 

in ro oo — n 04 

0 ) ■<* O — eo v 

c- t~ o a> cm ® 

• 

'OcsnaB 

-OBMOO 

— O co r~ co O 

— O <J) ® co O 

CD - - - 04 - 


-’-OOO'- 

--OOO- 

- - O O O - 

- - O O O - 



n 04 oo co in in 

04 O 00 04 o- O 

04 0- in ^ Tf TT 

M N O 0) O CM 

co in co ^ C7) 


in CD oo d- m ^ 

<7 S 00 CO — ID 

V - Tf CD CM O 

tr m co in 

M-OCMM 

o 

ip io oo in 0 ) *■ 

in in ^ in 04 oo 

— co in in 04 

in — O r- v t 

O’ in o* ^ o- in 

o 

^ o* co co ip co 

rr o- n o> co 

■n in 04 oo m © 

in ^ in co in 

P- — N O O V 


O 0) (7) 00 0) 0) 

O (7) C7) 00 (7) CD 

O (7) <7> h* 0) (7) 

O (7) (7) O* (7) 0) 

0) 0) GO 00 (7) O 


-ooooo 

-ooooo 

-ooooo 

-ooooo 

oooooo 

< 

o o o o o o 

oooooo 

oooooo 

oooooo 

oooooo 

Ui 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

z 

O d n in 

p- O co n in 

O co co in 

0* o co CO in 

O co co in 


04 n co 

04 co co 

04 (0 CO 

04 CO CO 

04 CO CO 


* 

04 

o 

0) 

r- 


o* 

in 

CO 

o 

® 

D 


IX) 


in 

o 



co 

CO 

04 

04 


6 

o 

6 

o 

o 



120 


TABLE V (Continued) 
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TABLE VII 


COMPUTER CODE LISTING FOR REDUCTION OF HOT-WIRE 

VOLTAGES 


COMPUTER PROGRAM TO C*LCULATF TURBULENCE 
QUANTITIES USING THE EXPERIMENTAL DAIA 
OBTAINED BV SIX ORIENTATION HOT-MIRE TECHNIQUE 


VERSION OF OCT. 1984 


PREPARED BV : 

SALIM I. JAN JUA 

MOOIFIED BY: 

GARY B FERRELL 

SCHOOL OF MECHANICAL AND AEROSPACE ENGINEERING 
OKLAHOMA STATE UNIVFRSITV 
STILLWATER OK. 74078 


01 MENS I M( 12).ER( 12 ) . AMECV! 12).VAR( 12) 

0 1 MENS I C (DUMO ( 6 ) . I IPDUMO ( 6 ) . VMOUMO ( 6 ) . VPDUMO ( 6 ) 

DIMENSION WMOUMO ( 6 ) . WPDUMO ( 6 ) . U VDUMO 16). UWDUMO C 6 ) . 

• VWDUMO ( 6 ) 

DIMENSION UMA ( 30 ) . VMA ( 30 ) . WMA ' 30 ) , UPA ( 30 ) , VPA ( 30 ) , 

• WPA ( 30 ) . UV$A< 30 ) . UWSA ( 30 ) . VWSA ( 30 ) . RADL ( 30 ) . KE ( 30 ) . AKF ( 30) 
REAL KE 

COMMON UBAR . VB AR , WBAR . UPRMS . VPRMS . WPRMS . 

•UVSS.UWSS. VWSS 
DATA 0IA.EITA/5.75.0 8/ 

DO B 7 IFILE-11. 16 
I OUT »IFILE*42 
REWINO IFILE 
JMAX "17 
IS-0 

READ( IFILE .* ) XDDIA, THETA. EMO 
RE AD( I F I LE . * ) A1.B1.C1 
READ( IFILE.*) A2.B2.C2 
READ( IFILE.*) A3.B3.C3 
WRITE! 10. 1412) 

1412 FORMAT! 1H1 ) 

WR1 tE( 1C ,1*1 1 • ) ''CIA 1 * *t T A 

1311 FORMAT!///. 4X. AXIAL POSITION. X/OIA • 4./. 

• 4 X , 'THETA • ' .F6. 1.//) 

WRITE! 10. 1111) 

1111 FORMAT!/. 4X. 'THE CALIBRATION CONSTANTS ARE:') 

WRITE! 10.*) AI.B'.CI 
WRITE! 10.*) A2.B2.C2 
WRITE! 10.*) A3.B3.C3 
X "XDD I A *0 1 A 
DO 77 JCOI INT • 1 . JMAX 
JC02-JC0UNT/2 

RF ADI IFILE . • WCEMU ) ,F»! I)). !• 1 .8) 


00000010 
00000020 
00000030 
00000040 
00000050 
00000060 
00000070 
00000080 
00000090 
OOOOOIOO 
000001 10 
00000120 
00000130 
00000140 
00000150 
00000160 
00000170 
00000180 
00000190 
00000200 
"00002 10 
00000220 
00000230 
00000240 
00000250 
00000260 
00000270 
00000280 
00000290 
00000300 
00000310 
00000320 
00000330 
00000340 
00000350 
00000360 
00000370 
00000380 
00000390 
00000400 
00000410 
00000420 
00000430 
00000440 
00000450 
00000455 
00000457 
''ww i c o 

00000470 

000004B0 

00000490 

00000500 

00000510 

00000520 

00000530 

00000540 

OOOOC550 

00000555 

f<'00>5«C 
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TABLE VII (Continued) 

I 

— ■ 


MM AL ( •- UCOUNT ) -0 . 500 
RDOI A-R/DI A 

IF ( ( JCD2*2 ) . CO JCOUNT ) WRITE ( 10. 1412) 

WRITE! 10. 1312 ) ROD I A 

1312 FORMAT!///. 4X. RADIAL POSITION. R/OIA • ’.FT «.//) 
WRITE! 10. 1112) 

1112 FORMAT!/. 7X. MEAN ANO R.M.S. VOLTAGES:'./) 

WRITE! 10. 1 100) (EM(I), 1-1.6) 

WRITE! 10. 1200) ( ER( I ) . 1-1.6) 

WRITE! 10. 112) 

1 100 FORMAT ( 6F9 . 4 ) 

1200 FORMAT (6F9 4) 

IS- IS^ 1 

RAOL! I S ) -RDOI A 

UM01-! -B1* SORT! B1*« 2-4 0*C 1 • ( A 1 - EMO* # 2 ) ) )/(2 0*C1 ) 
UM 1 • ( - B 1-SORT !B 1 • *2-4 . 0*C 1 • ( A 1 -FM( 1 ) • *2 ) ) )/( 2 0*C1 ) 
UMO -UMO 1 -UMO 1 
UM-UM1 -UM1 

DEU*B1/(4 . 0*EM( 1 )-UM1 )-C 1/( 2 0*EM( 1 ) ) 

UOEU-UM-DEU 
UPDUM-ER! 1 )/UDEU 
UMDUMO-UM/UMO 
UPDUMM-UPDUM-UMOUMO 
00 30 1-1,6 
EM2-EM( I )-EM( I ) 

ER2-ER! I )*ER( : ) 

O-SORT (B2**2-!4*C? , (A2-EM2) ) ) 

PHE • ( ! -B2*D)/(2*C2))**2 
DPHE — ( 2 • C M ( I )/C2)*( 1 - ( B2/0 ) ) 

D2PHE-I 1/EM( I ) )»0PHE*(B*B2*EM2)/0»»3 


LOCAL MEAN EFFECTIVE COOLING VELOCITY IS CALCULATED 


AMECV! I ) *PHF *0 5*D2PHE -ER2 


VARIANCE. VAR IS CALCULATED 


VAR ( I )-( (0PHE**2)*(ER2) *-((0 5-D2PHE •ER2)**2) 
AMECV! I-6)*AMECV< I ) 

VAR( I *6 )-V/.R( I ) 

WRITE! 10. 1 10) AMFCV( I ) . VAR! I ) 

110 r ORMAT( 7X . ' AMECV- ' . F7 . 4 , SX , ' VAR- ' . F7 4) 

30 CONTINUE 


MAIN CALLS THE SUBROUTINE CPYF TO CALCULATE 
THE PITCH ANO YAW FACTORS 


CALL CPYF( A 1 .B1.C1.A2.B2.C7. A3.B3.C3.PF, YF) 


PITCH FACTOR AND YAW FACTOR 


WRITE! 10.543) PF.YF 

543 FORMAT!/. 7X. 'PITCH FACT0R-'.F7 4.3X.'YAW FACTOR-'. F7 4) 
AL-PT-PF-YF-YF 
O-PF-PF-YF-YF 

WRITE! 10.444) UMOUMO . UPDUMM 

444 FORMAT!/. 7X. 'AXIAL MEAN VEL/INLET MAX VEL-'.FB 4.4X. 
•‘AXIAL TURB INTEN-'.FB 4) 

WRITE! 10.515) UMO 

515 FORMAT!/. 12X. MAX INLET VELOC I T Y • ' . F9 4 ) 

00 222 111-1,6 

II-I11-1 

N-6 

CALL STQTZ ( UDUMO . WMDUMO . VMOUMO . UPDUMO . WPOUMO . VPOUMO . 
•UVDUMO . UWDUMO . VWDUMO . N . I 1 1 ) 

C 


000005 IV 
00000580 
00000585 
00000590 
00000600 
00000610 
00000620 
00000630 
00000640 
00000650 
00000660 
00000670 
00000680 
00000690 
00000700 
00000710 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 
00000780 
00000790 
00000800 
00000810 
00000820 
00000830 
00000840 
00000850 
00000860 
00000870 
00000880 
00000890 
00000900 
00000910 
00000920 
00000930 
00000940 
00000950 
00000960 
00000970 
00000960 
00000990 
00001000 
OOOOIOIO 
00001020 
00001030 
00001040 
00001050 
00001060 
00001070 
00001080 
00001090 
00001 100 
OOOO’ 1 10 
00001120 
00001130 
00001140 
00001150 
00001160 
00001170 
00001180 
OOOOl 190 
00001200 
00001210 
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TABLE VII (Continued) 


MAIN CALLS the SUBROUTINE FMCV TO FIND THE 00001220 

THE MINIMUM COOLING VELOCITY AND THE TWO 00001230 

--ADJACENT ONES 00001240 

- 00001250 

CALL FMCV( AMECV N . I P , 1 0 . I R . I I ) 00001260 

ZP-AMECV(IP) 00001 2 7 0 

ZQ= AMECV( IQ) 00001280 

ZR=AMECV(IR) C0001290 

I F ( I 0 GT 6) 10=10-6 00001300 

I F ( I R . GT . 6 ) I R= ! R-6 00001310 

- - - 00001320 

-MAIN CALLS THE SUBROUTINE . SE ABC TO SET UP 00001330 

THE EOUATIONS FOR AO. BO. AND CO 00001340 

00001350 

CALL SEABC(ZP.ZQ.ZR.IP.AO.BO.CO) 00001360 

F = SQRT ( ( AO* * 2 ) * ( BO* * 2 )/ 3 ) 00001370 

I F ( CO . LT . F *0/AL ) GO TO 222 U0001380 

00001390 

VELOCITY FUNCTIONS F1.F2.AND F3 ARE CALCULATED 00001400 

- - - 00001410 

F1«S0RT( ( 1/AL)*( AO*F) ) 00001420 

IF( ( 1/AL)*( -AO*F ) . LT .0) GO TO 222 00001430 

F2=S0RT( ( 1/AL )*( -AO*F ) ) 00001440 

F3 = S0RT(C0-(0/AL )*F ) 00001450 

IF(F2 EQ.O) GO TO 222 00001460 

- - 00001470 

MAIN CALLS THE SUBROUTINE COABC TO CALCULATE 00001480 

THE FIRST AND SECOND DIFFERENTIALS OF AO. BO. 00001490 

AMD CO 00001500 

00001510 

CALL CDABC(DAP,DBP, OCP .D2AP.D2BP.D2CP.DAO, DBO .DCQ.D2AQ.D2B0, 00001520 

*D2CO.DAR,DBR.DCR.D2AR.D2BR.D2CR.ZP.ZO. ZR.IP) 00001530 

- MAIN CALCULATES THE FIRST AND SECOND 00001540 

DIFFERENTIALS OF THE VELOCITY FUNCTIONS 00001550 

F1.F2.AND F 3 WITH RESPECT TO THE 00001560 

SELECTED SET OF THE THREE COOLING VELO 00001570 

CITIES. 00001580 

- 00001590 

X 1 =F 1 *F 1 00001600 

X2 = X 1 *F 1 00001610 

X3=B0/(3*AL*AL) 00001620 

X4 = X 1/AL 00001630 

X5=(2*X2)-(2*F 1*A0/AL) 00001640 

X6=-(6»X1-2*A0/AL) 00001650 

Y1=F2*F2 00001660 

Y2*Y 1 *F2 00001670 

Y3=2.O*Y2*2.O*F2*A0/AL 00001680 

Y4=Y1/AL 00001690 

Y5=-(6*Y1*(2 0*A0/AL)) 00001700 

Z 1 *F3*F3 00001710 

Z2=Z 1 *F3 00001720 

Z3=2.0*Z2-2 0*C0*F3 00001730 

Z4=-(6 0* Z 1 - 2 0*C0) 00001740 

DF 1P=(X3*DBP*X4*DAP)/X5 00001750 

DF2P*(X3*DBP-Y4*DAP )/Y3 00001760 

DF3P* ( DCP* ( Z 1 -CO )♦ ( (0*0)/(AL*aL ) ) • ( A0*DAP+ ( BO*DBP ) /3 ) )/Z3 00001770 

DF 1Q«(X3*DBQ+X4*0AQ)/X5 00001780 

0F20*(X3*DB0-Y4*DA0)/Y3 00001790 

DF30* ( DCO* ( Z 1 -CO ) * ( (0*0)/(AL*AL ) ) * ( A0*DAQ* ( BO*DBO ) /3 ) )/Z3 00001800 

DF 1R.(X3*DBR+X4*DAR)/X5 00001810 

DF2R=(X3*DBR-Y4*DAR )/Y3 00001820 

DF3R-(DCR*(Z1-C0)*( (0*U)/(AL*AL ))*(A0*DAR»(B0*DBR)/3))/Z3 00001830 

02F IP" ( (X6*DF 1P*DF 1 P ) ♦ ( 2 .0*F 1/AL )*(DAP*DF 1P + '~AP*DF 1P)*(D2AP 00001840 

• *X 1/AL )♦ ( 1/(3*AL*AL) )*(DBP*DBP*B0*02BP) )/X5 00001350 

D2F2P«( ( Y5*DF2P*DF2P)-(2 0*F 2/ AL ) • ( DF 2P*DAP*OAP*OF2P ) - < Y 1 *D2 00001900 

♦AP/AL H f DBP*DBP+B0*02PP I/( 3 . 0*AL*AI ))/»3 0000*870 
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TABLE VII (Continued) 


D2F3P»(1 < E4*UFJF’*U*dKJ*2*»-J*1urjK •U(.K*uk.fur jrri uwr 'uu /’ui 
*-C0)*D2CP*( (0*0)/( AL*AL) )•( ( AO*02AP*0AP*DAP ) * ( DBP *DBP 
•*B0*D2BP)/3))/Z3 

02F 10= ( (X6*DF1Q*DF 1Q)*(2.0*F 1/AL ) • ( DAQ*DF 1Q*DAQ*0F 10)*(D2A0 
••X1/AL )♦( 1/(3.0*AL*AL) ) * ( DBQ*DBQ*B0*D2RQ ) ) / X5 
D2F20'( ( Y5*DF2Q*DF2Q)-(2.0*F2/AL)*(DF2Q*DAQ*DAQ*DF2Q)-( Y1*D 
*2 AO/ AL )♦ ( D80*0BQ*B0*D2BQ )/ (3.0*AL*AL ) )/Y3 
D2F3Q=( ( Z4*DF3Q*0F3Q)*2.0*F3*(DF3Q*DCQ*DCQ*DF3Q)-( DCQ’PCQ ( 

•♦ ( Z 1 -CO )*D2CQ*( (0*0)/(AL*AL ) )*< ( AO *02 AQ*DAQ»DAQ )♦ ( 
•DBQ*DBQ*B0*D2B0 )/3 ) )/Z3 

D2F 1R=( (X6*DF 1R*DF 1 R )♦ ( 2 . O* F 1/ AL ) • ( DAR»DF 1 R*DAR*DF 1 R ) + ( 02 AR 

* * X 1 / AL ) * ( 1/(3.0*AL*AL))*( 0BR*DBR*B0*D2BR ) )/X5 
D2F2R=((Y5*DF2R*DF2R)-(2.0*F2/AL)*(DF2R*DAR*DAR*DF2R)-(Y1*D 

•2AR/AL )♦ ( DBR*DBR*B0 *D2BR )/(3.0*AL*AL ) )/Y3 
D2F3R= ( (Z4*DF3R*DF3R)*2 . O* F 3* ( DF 3R*DCR*DCR»DF3R ) - ( OCR *DCR ) 
•*(Z1-C0)*D2CR*( (0*0)/(AL*AI ) )•( ( A0*D2 AR*DAR *0AR ) + ( 
»DBR*DBR*B0*D2BR), 3) )/Z3 

D2F 1PQ=( (XG‘DF 1P*DF 10 )♦ ( 2 . 0*F 1 / AL ) * ( 0AP*PF 1Q*0A0*DF IP ) * ( X 1 
•*D2APQ/AL)+( 1/(3 0‘AL‘AL) )‘(DBP‘DB0*B0‘D2BPQ) )/X5 
D2F 1QR = ( (XG*DF 10*0F 1 R ) ♦ ( 2 . 0*F I / AL ) * ( DAQ'OF 1R + DAR *DF 1 Q ) + ( X 1 

* *D2AQR/ AL )♦( 1/(3 0*AL*AL ) ) »(DBQ*DBR*B0*D2BQR ) )/X5 

D2F 1 PR= ( ( X6‘DF 1P*DF 1R)t(2.0*F 1 / AL ) • ( DAP ‘OF 1R+DAR*DF IP )♦( X 1 

* *02 APR/ AL ) * ( 1/(3.0*AL»AL) )*(DBP*DBR*B0*D2BPR ) )/X5 
D2F2PO=( ( Y5*DF2P*DF2Q)-(2.O*F2/AL)*(DF2P*DAQ*DAP*0F2Q)-(Y1 

••02APQ/AL )*( 1/(3. 0* AL*AL ) ) • ( OBP *DBQ*B0*U2BPQ ) ) / Y3 
02F20R=( ( Y5*DF20*DF2R)-(2 0*F2/AL )‘(DF20*0AR+DA0*DF2R)-( Y 1 
••02AQR/AL)M 1/(3.0*AL‘AL ) ) * ( 0BQ*0BR*B0*D2BQR ) ) / Y3 
D2F2PR=( (Y5‘DF2P‘DF2R)-(2.0‘F2/AL)‘(0F2P*0AR*0AP‘DF2R)-(Y1 
“D2APR/AL )♦( 1/(30* AL*AL ) ) • ( DBP *DBR*B0*D2BPR ) )/Y3 
D2F3P0- ( ( Z4*DF3P*DF3Q)+2 . 0* F 3 * ( DF3P ‘OCO+DCP *DF 30 > - ( DCP *DCQ ) 
**(Z1-C0)*D2CPQ*( (0*0)/(AL»AL ) ) * ( ( A0*D2 APO*OAP*DAO ) + ( 
*DBP*DBQ*B0»028PQ>/3) )/Z3 

02F30R* ( ( Z4*DF30*DF3R)*2 .0*F3*( 0F30*DCR*DCQ*0F3R )-(DCO*OOR) 
•♦(Z1-CU) *D2CQR* ( (0*0)/(AL*AL ) )*( (A0*D2AQR*DAQ*DAR)+( 
*DB0*DBR*B0*D2BQR)/3) )/Z3 

D2F 3PR= ( (Z4*DF3P*DF3R)*2 0*F 3* ( 0F3P *DCR+DCP*DF 3R ) - ( DCP ‘OCR ) 

* ♦ ( Z 1 -CO )*D2CPR+( (0*0)/(AL‘AL))*( ( A0*D2APR*DAP *DAR >* ( 

*DBP *DBR*BO *D2BPP )/3 ) )/23 


MAIN CALLS THE SUBROUTINE COVAR TO 
--CALCULATE THE COVARIANCE BETWEEN THE 
--SELECTED COOLING VELOCITIES. 


AKPQ=0.9*SQRT(VAR( IP ) *VAR( IP* 1 ) ) 

AKQR '0.9 ‘SORT ( VAR( IP*1 )*VAR( 1P+2)) 

AKPR=0 . 8 1 ‘EITA*SORT( VAR( IP )‘VAR( IP* 2) ) 

AKQP'AKPO 

AKRO'AKQR 

AKRP= AKER 


MAIN CALCULATES THE AX I AL . RADI AL . AND 
TANGENTIAL MEAN VELOCITIES 


UMEAN-F 1*0.5* (D2F 1P‘VAR( IP)*02F 1Q‘VAR( IP* 1 )*D2F 1R*VAR( IT+2 ) ) 
• + D2F 1PQ*AXPQ*D2F 1QR‘AK0R*D2F 1PR‘AKPR 
WME AN«F 2*0 .5*(D2F2P*VAR( I P ) + D2F 20* VAR( IP*1 ) +D2F 2R * VAR ( I P*2 ) ) 
*+D2F2P0‘AKP^*D2F20R‘AK0R*02F2PR‘AKPR 
VMEAN=F3*0.5*(D2F3P‘VAR( I P ) *D2F 30‘ VAR ( IP* 1 ) *02F3R * VAR( I P*2 ) ) 

•♦D2F3PQ‘AKPQ*D2F30R*AKQR*D2F3PR‘AKPR 

UP 1 =DF 1P‘0F 1P‘VAR( I P ) +DF 1Q*DF 10*VAR( IP* 1 ) + DF 1R*DF 1R‘VAR( IP 
** 2 ) 

UP2*DF 1P*DF IQ* AKPO*DF 1P»DF 1R*AKPR*DF 10* DF 1P*AK0P*DF 10*DF 1R»A 
•KOR+OF 1 R *0F IP* AKRP*DF 1R*0F IQ* AKRO 
UP3«0 . 5* ( D2F 1P*VAR(IP)*D2F10*VAR( IP*1 )*D2F 1 R* VAR( I P* 2 ) ) 
UP4=02F 1PQMKPQ + D2F 1QR*AK0R*02F 1PR* AKPR 
UP5*UP3*UP4 


00001890 
00001900 
00001310 
00001920 
00001930 
OOOC1940 
00001950 
00001960 
00001970 
00001980 
00001990 
00002000 
00002010 
00002020 
00002030 
00002040 
00002050 
00002060 
00002070 
00002080 
00002090 
00002 100 
000021 10 
00002120 
00002 1 30 
00002 1 40 
00002150 
00002 1 60 
00002 1 70 
00002 1 80 
00002 1 90 
00002200 
00002210 
00002220 
00002230 
00002240 
00002250 
00002260 
00007770 
00002280 
00002290 
00002300 
00002310 
00002320 
00002330 
00002340 
00002350 
00002360 
00002370 
00002380 
00002390 
00002400 
000024 10 
00002420 
00002430 
00002440 
00002450 
00002460 
00002470 
00002480 
00002490 
00002500 
00002510 
00002520 
00002530 
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TABLE VII (Continued) 


UPRMS2-UP1*UP2-UP5**2 

WP1-DF2P*0F2P*VAR< IP)*DF2Q*DF2Q*VAR( IP* 1 )*DF2R*OF2R*VAR( IP 

*♦2) 

WP2*DF2P*DF2Q*AKPQ*DF2P*DF2R*AKPR*DF2Q*DF2P*AKQP*DF20*0F2R*A 
*KQR*DF 2R *DF 2P* AKRP+DF2R *DF2Q* AKRO 
WP3'0.5*(D2F2P*VAR( I P ) +D2F20* VAR ( I P* 1 )*D2F2R*VAR( IP*2) ) 

WP4 'D2F 2PQ * AKPQ*D2r 2QR* AKQR*D2F 2PR* AKPR 

WP5*WP3*WP4 

WPRMS2*WP 1 *WP2 - WP5* *2 

VP 1»DF3P*DF3P*VAR( I P ) *DF 3Q*DF 3Q* VAR( IP* 1 )*DF3R«DF3R*VAR( IP 

• *2 1 

VP2*DF3P*DF3Q*AKPQ*CF3P*DF3R*AKPR*DF3Q*DF3P*AKQP*0F30*DF3R*A 
•KQR*DF3R*DF3P*AKRP+DF3R*0F3O»AKRQ 
VP3'0 5*(D2F3P*VAR( I P ) *D2F30* VAR ( IP* 1 ) + D2F3R*VAR( IP* 2 ) ) 

VP 4 'D2F 3PQ* AKPQ*D2F 3QR* AKQR + 02F 3PR • AKPR 

VP5'VP3*VP4 

VPRMS2=VP1+VP2-VP5**2 

UVI'DF 1P*DF3P*VAR( IP )*DF 10*DF 30* VAR ( I P* 1 )*DF 1R»DF3R*VAR( IP 

«*2) 

UV2'DF 1P*DF30*AKP0*0F 1P*DF3R*AKPR*0F 10*DF3P* AKQP*DF 1Q*DF3R*A 
•KQR + DF 1R*DF3P*AKRP + DF 1R'DF3Q*AKRQ 
UV3'0 . 5* ( D2F 1P*VAR( IP )*D2F 10* VAR ( IP* 1 )*D2F 1R* VAR( IP*2) ) 
UV4-D2F 1 PQ* AKPQ+U2F 1QR*AKQR*D2F 1PR»AKPR 

UV5'0.5*(02F3P*VAR( IP )*D2F30*VAR( IP*1 )*D2F3R»VAR( I P* 1 ) ) 

UV6*D2F3P0*AKP0*D2F3QR*AK0R*D2F3PR*AKPR 

UVPB'UV 1 *UV2 - ( ( UV3+UV4 )*(UV5*UVG) ) 

VW1'DF3P»0F2P*VAR( I P ) *DF30*DF20* VAR ( IP*1 )+DF3R*DF2R*VAR( IP 

**2 ) 

VW2'DF3P*DF2Q*AKPQ*DF3P*DF2R*AKPR*DF30*DF2P*AKQP*DF3Q*DF2R*A 
•K0R*DF3R*DF2P«AKRP*DF3R*0F20*AKR0 
VW3*0. 5*(02F3P*VAR( I P ) *D2F 30* VAR( I P* 1 )*02F3R‘VAP( IP* 2) ) 
VW4'D2F3PQ*AKPQ*02F30R*AKQR*D2F3PR*AXPR 

VW5'0. 5*( D2F2P*VAR( I P ) *D2F 20* VAR( I P* 1 ) *D2F 2R * VAR( I P* 1 ) ) 

VWG -D2F2PQ* AKPQ*D2F 2QR* AK0R + D2F 2PR*AKPR 
VWPB= VW 1 * VW2 - ( ( VW3*VW4 )• ( VW5+VW6 ) ) 

UW 1 =DF IP *0F 2P * VAR ( IP)*OF 1Q*DF 20* VAR( I P* 1 )*0F 1R'DF2R*VAR( IP 

•*2 ) 

UW2'0F 1 P*OF 20* AKPO*OF 1P*DF2R*AKPR*DF 10*0F2P* AKQP*0F 1Q*DF2R*A 
•KQR+DF 1R*DF2P*AKRP+DF 1R*DF20*AKR0 
UW3'0 5* ( D2F 1P*VAR( I P ) *D2F 10* VAR ( I P* 1 ) + 02F 1R*VAP( IP*? ) ) 
UW4=D2F 1PQ»AKP0*02F 1QR*AKQR*02F 1PR*AKPR 
UW5'0 5*(D2F2P*VAR( IP )*0?F20 *VAR( I P+ 1 ) *D?F 2R* VARf I P+ 1 ) ) 
UW6-D2F2P0*AKP0*D2F20R*AK0R*P2F2PR*AKPR 
UWPB'UW 1 + UW2 - ( ( UW3+UW4 ) • ( UW5+UW6 ) ) 

UDUMO ( I 1 1 ) 'UME AN/tIMO 
WMDUMO ( I 1 1 ) = WME AN/UMO 
VMDUMOl I 1 1 ) =VMf AN/UMO 

I F (UPRMS2 . GT . 0 . O) UPOUMOI I 1 1 ) 'SORT ( UPRMS2 ) /UMO 
I F ( WPRMS2 . GT 0 0) WPDUMO( I 1 1 ) 'SORT ( WPRMS2 ) /UMO 
IF(VPRMS2 GT .0.0) VPDUMO ( 1 1 1 ) 'SORT ( VPRMS2 )/UMO 
UVDUMOl I 1 1 ) 'UVPB/UMO* • 2 
VWDUMO ( I 1 1 ) 'VWPB/UMO* *2 
UWDUMOI I 1 1 )*UWPB/UM0**2 


MAIN CALLS THE SUBROUTINE AVRG TO COMPUTE AN 

ENSEMBLE AVERAGE OF THE TIME-MEAN AND 

TURBULENCE QUANTITIES 


222 CONTINUE 

CALL AVRG( UDUMO . VMDUMO . WWOUMO . UPDUMO . VPDUMO . WPDUMO . 
•UVDUMO . UWDUMO . VWDUMO . N ) 

UMA (IS) *UBAR 
VMA ( IS) ■ VBAR 
WMA (IS) »WBAR 
UPA ( IS) 'UPRMS 
VF A ( 1 S ) *VPRMS 


00002540 
00002550 
00002560 
00002570 
00002580 
00002590 
00002600 
00002610 
00002620 
00002630 
00002640 
00002650 
00002660 
00002670 
00002680 
00002690 
00002700 
00002710 
00002720 
00002730 
00002740 
00002750 
00002760 
00002770 
00002780 
00002790 
00002800 
00002810 
00002820 
00002830 
00002840 
00002850 
00002860 
00002870 
00002880 
00002890 
00002900 
00002910 
00002920 
00002930 
00002940 
00002950 
00002960 
00002970 
00002980 
00002990 
00003000 
00003010 
00003020 
0000303 O 
00003040 
00003050 
00003060 
00003070 
00003080 
00003090 
00003 <00 
00003110 
00003120 
00003130 
00003140 
00003150 
00003160 
00003170 
00003180 
00003190 
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TABLE VII (Continued) 



WPA( IS) -WPRMS 

UVSA(IS)-UVSS 
UWSA( IS)*UWSS 
VWSA(IS)*VWSS 

AK* ( ( UP RMS * *2 ) ♦ ( VPRMS • *2 ) ♦ ( WPRMS * *2 ) )/2 
KE ( IS ) * AK 
AKE (IS ) *UPOUMM 
1 12 FORMAT ( ' ' ) 

WR I TE ( 10.112) 

WRITE( 10. 1000) UBAR.VFIAR.WBAR 

WR I TE ( 10.2000) UPRMS. VPRMS. WPRMS. AK 

WRITE( 10. 3000) UVSS.UWSS. VWSS 

lOOO FORMAT (/. ' UBAR- ' . F7 4 . 7X . ' VBAR* ' . F7 . 4 . 7X . ' WBAR= ' , F7 . 4 > 

2000 FORMAT ( / , ' UPRMS* ' . F7 4 . 7X . ' VPRMS = ' . F 7 4 . 7X . 

•'WPRMS*' .F7.4.7X, *K-'.F7.4) 

3000 FORMAT (/. ' UVSS* ' . F 7 . 4 . 7X . ' UWSS*' .F7.4.7X. ' VWSS* ' . F7 . 4 ) 
77 CONTINUE 
999 CONTINUE 


DO 127 I * 1 . JMAX 

WRITE( IOUT. 129) RADL ( I ) . UMA ( I ) . VMA ( I ) . WMA ( I ) , UP A ( I ) , 
1 VPA ( I ) .WPA( I ) ,UVSA( I ) , UWSA ( 1 ) . VWSA ( I ) . KE ( I ) . AKF ( I ) 

127 CONTINUE 

129 F0RMAT(4F9.5./.4F95./.4F9 5) 

87 CONTINUE 
STOP 
END 


THIS SUBROUTINE SETS TURBULENT QUANTITIES TO 
ZERO AT THE BEGINING OF EACH ITERATION 


SUBROUT I NE STQT Z ( UOUMO . WMDUMO . VMDUMO . UPDUMO . WPDUMO . 

* VPDUMO . UVOUMO . UWOUMO . VWOUMO . N . I ) 

DIMENSION UDUMO ( 6 ) , WMDUMO ( 6 ) . VMDUMO ( 6 ) . UPDUMO ( 6 ) . WPDUMO ( 6 ) 
DIMENSION VPDUMO (6 ) ,UVDUMO(G ) .UWDUMO(G) . VWOUMO (6 ) 

UDUMO ( I ) =0 . O 
WMDUMO ( I ) *0 0 
VMDUMO ( I ) *0 . 0 
UPDUMO( I ) =0 . 0 
WPDUMO( I ) *0 . 0 
VPDUMO ( I ) *0 0 
UVDUMOI I ) *0 . 0 
UWOUMO ( I ) *0 . 0 
VWDUMO ( I ) *0 . 0 
RETURN 
END 


THIS SUBROUTINE FINDS THE MINIMUM MEAN EFFECTIVE 
COOLING VELOCITY AND THE TWO ADJACENT TO IT 


SUBROUTINE FMCV(CV.N. IX.IY.IZ.II) 


WWJiUVJ 

00003210 
00003220 
00003230 
00003240 
00003250 
00003260 
00003270 
00003280 
00003290 
00003300 
00003310 
00003320 
00003330 
00003340 
00003350 
00003360 
00003370 
00003380 
00003390 
00003400 
000034 10 
00003420 
00003430 
00003440 
00003450 
00003460 
00003470 
00003400 
00003490 
00003500 
00003510 
00003520 
00003530 
00003540 
00003550 
00003560 
00003570 
00003580 
00003590 
00003600 
00003610 
00003620 
00003630 
00003640 
00003650 
00003660 
00003670 
00003680 
00003690 
00003700 
00003710 
00003720 
00003730 
00003740 
00003750 
00003760 
00003770 
00003780 
00003790 
00003800 
000038 10 
00003820 
00003830 
00003840 
00003850 
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TABLE VII (Continued) 



DIMENSION CV( 50) 


00003860 


IF ( CV( 2 ) . LT ,CV( 1)) GO TO 

20 

00003870 


I F ( CV( 3 ) LT . CV( 1 ) ) GO TO 

30 

00003880 


I F ( CV( 4 ) LT . CV( 1 ) ) GO TO 

40 

00003890 


I F ( CV( 5 ) . LT . CV( 1)) GO TO 

50 

00003900 


IF(CV(6) LT CV( 1 ) ) GO TO 

60 

00003910 


I X = 6 


00003920 


I V- 1 


00003930 


12*2 


00003940 


GO TO 100 


00003950 

20 

I F ( CV( 3 ) LT CV(2) ) GO TO 

30 

00003960 


IF(CV(4).LT.CV(2)) GO TO 

40 

00003970 


IF(CV(5) LT CV ( 2 ) ) GO TO 

50 

00003980 


IF(CV(6) .LT . CVI 2 ) ) GO TO 

60 

00003990 


I X - 1 


00004000 


I Y = 2 


00004010 


I Z * 3 


00004020 


GC TO 100 


00004030 

30 

IF(CV(4) LT.CV(3) ) GO TO 

40 

00004040 


IF ( CV( 5 ) LT CV( 3 ) ) GO TO 

50 

00004050 


IF(CV(6).LT.CV(3)) GO 10 

60 

00004060 


I X * 2 


00004070 


I Y * 3 


00004080 


IZ = 4 


00004090 


GO TO 100 


00004100 

40 

IF(CV(5).LT.CV(4)) GO TO 

50 

00004110 


I F ( CV( 6 ) LT CV( 4 ) ) GO TO 

60 

00004 1 20 


I X» 3 


00004 1 30 


I Y *4 


00004140 


I Z *5 


00004 150 


GO TO 100 


00004160 

50 

I F ( CV( 6 ) LT.CV(6) ) GO TO 

60 

00004 1 70 


IX'4 


00004180 


I Y *5 


00004190 


I Z*6 


00004200 


GO TO 100 


00004210 

60 

I X *5 


00004220 


I Y » 6 


00004230 


IZ* 1 


00004240 

100 

IXMX+II 


00004250 


I F ( I X GT 6) I X * I X -6 


00004260 


I F ( I Y . GT . 6 ) IY-IY-6 


00004270 


IFIIZ.GT 6) I Z= I Z-6 


00004280 


IY»IX+ 1 


00004290 


I Z* I X*2 


00004300 


RETURN 


00004310 


END 


00004320 




00004330 

* * * * 



00004340 




00004350 


THIS SUBROUTINE CALCULATES THE PITCH AND YAW 

00004360 


FACTORS USING THE THRE E -D I RECT IONAL CALIBRATION 

00004370 


CONSTANTS . 


00004380 




00004390 




00004400 




00004410 


SUBROUTINE CPYF ( A 1 , B 1 . C 1 . 

A2.B2.C2. A3.B3.C3.PF ,YF » 

00004420 


E = 3 . 0 


00004430 

10 

W1*B3**2-4 0*C3*(A3-E**2) 


00004440 


IF(WI.LT.O.O) GO TO 20 


00004450 


E*E*0 05 


00004460 


GO TO 10 


00004470 

20 

E-E-0.05 


00004480 


W1»( -B3*S0RT(B3‘»2-4 0*C3 

• ( A3 -E * *2 ) ) )/(2 .0*C3) 

00004490 


W-W1»W1 


00004500 


V1-( -B2 + S0RT ( B2* • 2 - 4 . 0*C2 

•(AJ-F^jm/U.O'CJI 

00004510 
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TABLE VII (Continued) 




V»V1*V1 

WWW’ 



U 1 - ( -B1*SQRT(B1*»2-4 0*C1*(A1-E**2) ) )/(2 0*C1 ) 

00004530 



U=U1*U1 

00004540 



PF*V/U 

00004550 



YF*V/W 

00004560 



RETURN 

00004570 



END 

00004580 

c 



00004590 

C ' 



00004600 

c 



00004610 

c 


THIS SUBROUTINE SETS EOUATIONS FOR AO. BO. AND CO 

00004620 

c 


DEPENDING UPON THE SET OF THE THREE COOLING 

00004630 

c 


VELOCITIES CHOSEN 

00004640 

c 



00004650 

C ' 



00004660 

c 



00004670 



SUBROUTINE SE ABC( A 1 . A2.A3.K.X.Y.Z) 

00004680 



I F ( K EO 1 ) GO TO 15 

00004690 



I F ( K . EQ . 2 ) GO TO 25 

00004700 



I F ( K . EQ . 3 ) GO TO 35 

000047 10 



IF(K EO 4) GO TO 45 

00004720 



IF(K.EO 5) GO TO 55 

00004730 



I F ( K . EO . 6 ) GO TO G5 

00004740 


15 

X* A2 * * 2 -A3* *2 

00004750 



Y « - 2 0*A1**2*3.0*A2**2-A3**2 

00004760 



Z«A1»»2-A2**2*A3**2 

00004770 



GO TO 105 

00004780 


25 

X* A 1 • • 2- A2 * *2 

00004790 



Y=-(A1**2)*3.0*A2**2-2.0*A3**2 

00004800 



Z«A1**2-A2**2*A3**2 

00004810 



GO TO 105 

00004820 


35 

X*A1**2-2.0*A2**2*A3**2 

00004830 



Y*A 1 • *2- A3* *2 

00004840 



Z*A1**2-A2**2*A3**2 

00004850 



GO TO 105 

00004860 


45 

X«-(A2**2)+A3**2 

00004870 



Y = - 2 0*A1**2*3.0*A2**2-A3**2 

00004880 



Z = A 1 • *2 - A2 * *2* A3* *2 

00004890 



GO TO 105 

00004900 


55 

X--(A1**2)+A2**2 

00004910 



Y*-(A1**2)*3.0*A2**2-2.0*A3**2 

00004920 



Z»A1**2-A2**2*A3**2 

00004930 



GO TO 105 

00004940 


G5 

X*-(A1**2)*2.0*A2**2-A3**2 

00004950 



Y * - ( A 1 • • 2 )*A3* *2 

00004960 



Z - A 1 **2-A2**2+A3» *2 

00004970 


105 

RETURN 

00004980 



END 

00004990 

c 



00005000 

c # 



000050 1 0 

c 



00005020 

c 


THIS SUBROUTINE CALCULATES THE FIRST AND SECOND 

00005030 

c 


DIFFERENTIALS OF THE FUNCTIONS AO, BO. AND CO WITH 

00005040 

c 


RESPECT TO THE THREE CHOSEN MEAN EFFECTIVE COOLING 

00005050 

c 


VELOCITIES 

00005060 

c 



00005070 

c • 



00005000 

c 



00005090 



SUBROUTINE CDAPC ( A 1 . B 1 . C 1 . A2 1 .B2 1 . C2 1 . A2 , B2 . C2 . A22 . B22 . C22 , 

00005100 



*A3.B3.C3, A23.B23.C23.X. Y. Z.K) 

00005110 



I F ( K . EQ 1 ) GO TO 16 

00005120 



IF|K EQ 2) GO TO 26 

00005 1 30 



IF(K EQ 3) GO TO 36 

00005140 



I F ( K . EQ . 4 ) GO TO 46 

00005150 



IF(K EO 5) GO TO 56 

00005160 



IF(K EO 6) GO TO 66 

00005170 


TABLE VII (Continued) 


A1-0.0 

00005 mo 

b i •-4*x 

00005190 

C1*2*\ 

00005200 

A21-0.0 

00005210 

B2 1*-4 .0 

00005220 

C2 1*2.0 

00005230 

A2*2 0*Y 

00005240 

B2*6 . 0* Y 

00005250 

C2* -2 . 0* Y 

00005260 

A22*2 .0 

00005270 

B22-6.0 

00005280 

C22*-2 .0 

00005290 

A3*-2.0*Z 

00005300 

B3*-2 .0*2 

00005310 

C3«2 .0*2 

00005320 

A23*-2 0 

00005330 

B23* - 2 . 0 

00005340 

C23-2 0 

00005350 

GO TO 106 

00005360 

A 1 *2 . 0*X 

00005370 

B 1 *-2 0*X 

00005380 

C 1 *2 . 0*X 

00005390 

A21*2 .0 

00005400 

B2 1 * -2 . 0 

00005410 

C2 1 *2 . 0 

00005420 

A2*-2 . 0*Y 

00005430 

B2*6 . 0* Y 

00005440 

C2*-2 .0*Y 

00005450 

A22--2 .0 

00005460 

B22*6 .0 

00005470 

C22--2 .0 

00005480 

A3-0 0 

00005490 

B3»-4 .0*2 

00005500 

C3*2 0*2 

00005510 

A23-0 

00005520 

B23--4 .0 

00005530 

C23-2 .0 

00005540 

GO TO 106 

00005550 

A1*2.0*X 

00005560 

B 1 *2 0*X 

00005570 

C1*2 0*X 

00005580 

A21-2.0 

00005590 

B21*2 .0 

00005600 

021*2.0 

00005610 

A2 ■ -4 . 0* Y 

00005620 

B2-0.0 

00005630 

C2*-2.0*Y 

00005640 

A22*-4 . 0 

00005650 

B22-0.0 

00005660 

C22«-2 .0 

00005670 

A3* 2 . 0* Z 

00005680 

B3--2 ,0*Z 

00005690 

03*2.0*2 

00005700 

A23*2 . 0 

00005710 

B23--2 .0 

00005720 

023*2.0 

00005730 

GO TO 106 

00005740 

A 1 *0 . 0 

00005750 

B 1 * -4 . 0*X 

00005760 

01*2 0*X 

00005770 

A2 1 *0 . 0 

00005780 

B21--4 .0 

00005790 

021*2.0 

00005800 

A2*-2 0*V 

00005810 

B2-fl 0*V 

OOUC88JO 

02* -2 . 0* Y 

00005830 
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TABLE VII (Continued) 


A22--2.0 


00003840 

622-60 


00005850 

C22--2 0 


00005860 

A3-2 0*Z 


00005870 

B3--2 0*Z 


00005880 

C3-2 .0 *Z 


00005890 

A23-2 . C 


00005900 

623- -2 . 0 


00005910 

C23-2.0 


00005920 

GO TO 106 


00005930 

A 1 * -2 . 0*X 


00005940 

B 1 - -2 ,0*X 


00005950 

C 1*2 .0*X 


00005960 

A21--2.0 


00005970 

B21--2.0 


00005980 

C21-2.0 


00005990 

A2*2 . 0* Y 


00006000 

62*6 .0* Y 


00006010 

C2--2 0*Y 


00006020 

A22-2.0 


00006030 

622*6.0 


00006040 

C22--2.0 


00006050 

A3-0.0 


00006060 

B3--4 0*Z 


00006070 

C3*2.0*Z 


00006080 

A23*0 . 0 


00006090 

B23=-4 .0 


00006100 

C23*2 .0 


00006110 

GO TO 106 


00006 1 20 

A 1 * -2 . 0* X 


00006130 

B1*-2.0*X 


00006140 

C 1-2 0*X 


00006150 

A21--2 .0 


00006160 

B21--2 .O 


00006170 

C2 i * 2 . 0 


'00006180 

A2*4 . 0* Y 


00006190 

B2-0.0 


00006200 

C2 *- 2 . 0* Y 


00006210 

A22*4 0 


00006220 

B22-0.0 


00006230 

C22--2 O 


00006240 

A3*-2 .0*Z 


00006250 

B3-2 0*Z 


00006260 

C3*2 . 0* Z 


00006270 

A23--2.0 


00006280 

B23-2 .0 


00006290 

C23-2 .0 


00006300 

RETURN 


00006310 

END 


00006320 



00006330 



00006340 



00006350 



00006360 



00006 370 



00006380 

THIS SUBROUTINE CALCULATES 

THE COVARIANCES BETWEEN THE 

00006390 

VELOCITY FLUCTUATIONS USING 

A METHOD SUGGESTED BY KING 

00006400 



00006410 



00006420 



00006430 



00006440 



00006450 



00006460 



00006470 

SUBROUTINE COVAR ( CV . V , N . I P . 

7P . ZD. ZR . AKPQ . AKPR . AKOP . AKOP . AKRP 

00006480 

.AKRQ.FITA) 


00006490 
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TABLE VII (Continued) 


DIMENSION CV( 50 ) , V( 50 ) 

EITA-0.8 
DO 15 I*1,S 
CV( I+6)-CV( I ) 

V( I+6)»V( I ) 

15 CONTINUE 

IF(V( IP) .LE .0.002 ) GO TO 108 

ZET A 1 *SQRT ( 7P* *2 -2 .0* ZQ* *2+2 . 0*ZR* *2 ) 

ZETA3=SQRT (2.0*ZP**2-2.0*ZQ**2+ZR**2) 

PI 1*CV( IP+3)-ZETA 1-0.5* ( ( 1/CV( I P + 3 ) - ZP ♦ « 2/CV( I P+ 3 ) * *3 ) * V( I P ) 
•-(4 0*ZQ**2/CV( I P+3) **3+2 0/CV( I P+3 ) )*V( IP* 1 )♦< -4 0*ZR*»2 
*/CV( I P + 3)* *3+2 0/CV( I P + 3) )«V( I P + 2) ) 

PI3»CV( IP*5)-ZET A3-0. 5*( (2 0/CV( I P + 5) -4 0*ZP**2/CV( IP*5)**3 
•)*V( IP )♦( -2 . 0/CV( I P+5 ) -4 . 0*Z0* *2/CV( IP+5 )• *3 ) »V( IP* 1 )♦< 1/CV( 
•IP+5)-ZR**2/CV( IP + 5)**3)*V( IP + 2 ) ) 

A1«-2.0*ZP**2*EITA/V( IP+1 ) 

B1*6 0*ZP*ZO-( ZP*EITA/(Z0*V( IP* 1 ) ) )*(PI 1*CV( I P + 3 ) * • 3-P 1 3 *CV 
•( IP+5 ) • *3 ) 

Cl'PI 1*CV( IP + 3)**3-2 . 0*PI3*CV( I P + 5)* *3 
IF(B1**2-4.0*A1*C1 LT .0) GO TO 57 
AKP01*( -B1+S0RT(B1**2-4 ,0*A1*C1 ) )/(2 0*A1 ) 

AKPQ2«( -B1 -SORT (01* *2 -4 O* A 1 *C 1 ) )/( 2 . 0* A 1 ) 

RPQ1»AKPQ1/SQRT( V( IP)*V( IP* 1 ) ) 

RPQ2* AKPQ2/SQRT ( v ( IP)*V( IP+1 ) ) 

I F ( ABS( RPO 1 ) . GT . 1 ) GO TO 17 
GO TO 27 

17 I F ( ABS ( RPQ2 ) . GT . 1 ) GO TO 37 
AKPQ* AKPQ2 
27 AKPQ'AKPQI 
GO TO 47 

37 AKP0*O.9*S0RT(V(IP)*V( IP+1 )) 

47 AKQR" ( 2 . 0*CV( IP )*CV( IP+1 )*KPO + PI 1*CV( I P + 3 ) ♦ *3-P I 3*CV( I P + 5 ) 

• * *3 )/( 2 . 0*CV( IP* 1 )*CV( IP + 2) ) 

RQR= AKQR* SORT (V( I P+ 1 ) • V( IP+2 ) ) 

I F ( ABS(ROR) GT . 1 ) AKQR *0 . 9*SQRT ( V( IP* 1 ) • V( * P *2 ) ) 
AKPR*EITA*AKPQ«AKQR/V( IP* 1 ) 

GO TO 107 

57 AKPQ*0 . 9* SORT (V(IP)*V(IP+1)) 

AK0R'0.9*SQRT(V(IP+1)*V( IP+2) ) 

AKPR’EIT A * AKPO* AKQR/V( IP+ 1 ) 

GO TO 107 
108 AKPQ-0 0 
AKQR*0 ■ 0 
AKPR*0 . 0 
107 AKQP*AKPQ 
AKRO+AKQR 
AKRP * AKPR 
RETURN 
END 


THIS SUBROUTINE CALCULATES THE ENSEMBLED AVERAGE 
OF THE TIME-MEAN AND TURBULENCE QUANTITIES 


SUBROUT I NE AVRG( U . V . W . UPR . VPR , WPR . UV . UW . VW . N ) 

DIMENSION U(6).V(6) . W( € ) . UPR ( 6 ) . VPR( 6 ) . WPR ( 6 ) . UV( 6 ) ,UW(6). 
• VW( 6 ) , IN( 9 ) 

COMMON UB AR . VBAR . WBAR . UPRMS . VPRMS . WPRMS , UVSS . UWSS . VWSS 

UBAR *0 

VBAR+O 

WBAR+O 

UPRMS-0 

VPRMS-0 

WPRMS-0 

UVSS-0 

UWSS-0 

VW5S-0 


00006500 
00006510 
00006520 
00006530 
00006540 
00006550 
00006560 
00006570 
00006580 
00006590 
00006600 
000066 10 
00006620 
00006630 
00006640 
00006650 
00006660 
00006670 
00006680 
00006690 
00006700 
00006710 
00006720 
00006730 
00006740 
00006750 
00006760 
00006770 
00006780 
00006790 
00006800 
000068 10 
00006820 
00006830 
00006840 
00006850 
00006860 
00006870 
00006880 
00006890 
00006900 
00006910 
00006920 
00006930 
00006940 
00006950 
00006960 
00006970 
00006980 
00006990 
00007000 
00007010 
00007020 
00007030 
00007040 
00007050 
00007060 
00007070 
00007080 
00007090 
00007 100 
00007110 
00007120 
00007 1 30 
00007140 
00007 150 
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TABLE VII (Continued) 



DO 10 1*1,9 

00007160 


IN< I)-0 

00007170 

10 

CONTINUE 

00007180 


DO 20 U* 1 .6 

00007190 


IF(U( J) LE 0) GO TO 30 

00007200 


UBAR*UBAR*U( J ) 

00007210 


IN( 1 )*IN( 1 )♦ 1 

00007220 

30 

IF(V( J) LE 0) GO TO 40 

00007230 


VBAR* VBAR* V( J ) 

00007240 


IN( 2 )■ IN( 2 )♦ 1 

00007250 

40 

IF(W( J) LE .0) GO TO 50 

00007260 


WBAR*WBAR*W( J) 

00007270 


IN( 3) ■ IN( 3 )♦ 1 

00007280 

50 

IF(UPR( J) LE .0) GO TO 60 

00007290 


UPRMS*UPRMS*UPR( J ) 

00007300 


INI 4 )■ IN( 4 )♦ 1 

00007310 

60 

I F ( VPR ( J ) LE 0) GO TO 70 

00007320 


VPRMS*VPRMS+VPR( J ) 

00007330 


IN(5)*IN(5)+1 

00007340 

70 

I F ( WPR( J ) LE 0) GO TO 80 

00007350 


WPRMS*WPRMS + WPR( J) 

00007360 


1N(6)*IN(6)*1 

00007370 

80 

IF(UV( J) . LE .0) GO TO 90 

00007380 


UVSS’UVSS*UV( J ) 

00007390 


IN(7 )» INC 7 )♦ 1 

00007400 

90 

I F ( UW( J ) . LE 0 ) GO TO 100 

000074 10 


UWSS*UWSS*UW( J) 

00007420 


IN( 8 ) ■ INI 8 )♦ 1 

00007430 

100 

I F ( VW( d ) . LE 0) GO TO 20 

00007440 


VWSS*VWSS*VW( J ) 

00007450 


IN( 9 ) * INI 9 )♦ 1 

00007460 

20 

CONTINUE 

00007470 


DO 25 1*1,9 

00007480 

25 

IFIIN(I).EO.O) I N< I ) — lOOOO 

00007490 


UBAR=UBAR/IN( 1 ) 

00007500 


VBAR* VBAR/ I N ( 2 ) 

00007510 


WBAR*WBAR/INI 3 ) 

00007520 


UPRMS*UPRMS/IN(4 ) 

00007530 


VPRM5*VPRMS/IN(5) 

00007540 


WPRMS*WPRMS/IN(6) 

00007550 
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Figure 5. Dilution Jet. 
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Figure 7. The Three-Dimensional Hot-Wire Calibration. 
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Figure 8. Pitch and Yaw Factors as a Function of Hot-Wire f 
Effective Voltage. 
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Figure 12. 


Data Acquisition 


and Probe Drive System. 
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Figure 13. Bubble Generator and Injection Setup. 
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Figure 14. Smoke Generator and Injection Setup. 
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Spark-Gap Equipment Schematic. 
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Figure 17. 


Smoke Flow Visualization for Jet-to-Crossf low 
Velocity Ratio P = ?.G, 4.0, 6.0. 
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Figure IB. 9'jbble How Visualization f ur Jet-to-Crossf low 
Velocity Ratio R = ?.0, Swirl Vane 
: = 0, 45, 70. 
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Figure 19. Bubble Flow Visualization for Jet-to-Crossf low 
Velocity Ratio R = 4.0, Swirl Vane Angle 
4> - 0, 45, 70. 
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Figure 20. Bubble Flow Visualization for Jet-to-Crossf low 
Velocity Ratio R = 6.0, Swirl Vane Angle 
4> = 0, 45, 70. 
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Figure 21. Spark-Gap Flow Visualization for Jet-to-Crossf low 
Velocity Ratio = 2.0, Swi>-1 Vane Angle 
* = 0, 45, 70. 
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Figure 22. Spark-Gap Flow Visualization for Jet-to-Crossf low 
Velocity Ratio = 4.0, Swirl Vane 
Angle 0 = 0, 45, 70. 
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Figure 23. Spark-Gap Flow Visualization for Jet-to-Crossf low 
Velocity Ratio = 6.0, Swirl Vane 
Angle = 0, 45, 70. 
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Figure 24. Time-Mean and Turbulent Flowfield, R = 2.0, Traverse 
Angle P = 270 Degrees . 
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Angle 0 = 330 Degrees. 
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Figure 28. Time-Mean and Turbulent Flowfield, R = 2.0, Traverse 
Angle 0 = 30 Degrees. 
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Figure 28. (Continued) 
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Figure 29. (Continued) 
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Figure 30. Time-Mean and Turbulent Flowfield, R = 
Angle G = 270 Degrees. 
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Figure 30. (Continued) 
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Figure 31. Time-Mean and Turbulent Flowfield, R = 4.0, Traverse 
Angle 0 = 300 Degrees. 
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Figure 31. (Continued) 
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Figure 32. (Continued) 
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Figure 34. (Continued) 
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Figure 35. Time-Mean and Turbulent Flowfield, R = 4.0, Traverse 
Angle 0 = 60 Degrees. 
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Figure 37. Time-Mean and Turbulent Flowfield, R = 6.0, Traverse 
Angle 6 = 300 Degrees. 
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Figure 39. (Continued) 
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Figure 40. Time-Mean and Turbulent Flowfield, P. = 6.0, Traverse 
Angle 0 = 30 Degrees. 
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Figure 41. Time-Mean and Turbulent Flowfield, R = 6.0, Traverse 
Angle 6 = 60 Degrees. 
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